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ABSTRACT 


This thesis investigates a means to mitigate co-channel interference from fourth generation 
cellular signals in order to support passive synthetic aperture radar (SAR) imaging using 
those same signals. Synthetic aperture radar is a staple of the military and intelligence 
communities, but the active transmission required for such images informs the target of the 
imaging process. Use of passive signals, such as the orthogonal frequency-division multi¬ 
plexing (OFDM) signals of fourth generation cellular systems, is an attractive option, but 
co-channel interference mitigation is required. A method to separate the transmitted sig¬ 
nals that leverages the estimated signal delays between multiple transmitters and receivers 
is examined for narrowband and wideband signals. Multiple methods of assessing recovery 
performance are proposed. The impact of noise is considered, as is the impact of collection 
geometry on recovery performance, and recovery of OFDM signals is evaluated. Signal 
interpolation is a critical element of the proposed recovery process, and two methods are 
compared for accuracy and speed of computation. 
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Executive Summary 


Introduction 

Imaging of the earth by synthetic aperture radar (SAR) provides a number of advantages 
over the more traditional electro-optical methods, as SAR operates without a requirement 
for daylight illumination of the target area and can image through clouds, providing visi¬ 
bility to targets that would not otherwise be possible to image [1], Traditional monostatic 
SAR systems require illumination of the target area with a transmitter, which, in the case 
of military or intelligence targets, can serve to alert the target area to the imaging process. 
This affords the opportunity to shelter equipment or obscure the target area with an active 
radar jammer, negating some of the benefits of the SAR imaging approach [2]. 

Passive radar systems exploit signals that already exist in the target area, mitigating the 
warning provided by the monostatic SAR system [3]. The proliferation of fourth genera¬ 
tion cellular communications networks coupled with the performance of their orthogonal 
frequency-division multiplexing (OFDM) waveform makes those transmitters an appealing 
source to exploit for passive SAR imaging. 

Unfortunately, that same proliferation poses difficulty to the SAR imaging process from 
the satellite-based receivers in the form of co-channel interference. Cellular networks rely 
upon spatial frequency reuse in order to make the best possible use of the limited spectrum 
available. By dividing a larger geographic area into smaller cells and assigning each cell a 
subset of the available frequencies, it is possible to ensure that nearby cells do not reuse a 
particular frequency. Satellite-based receivers view many of these transmitters simultane¬ 
ously, and a means other than line-of-sight is required to isolate the individual transmitter 
signals. 


Recovery Process 

Consider a small system, consisting of two transmitters (transmitting the signals ■s'i(t) and 
S 2 (t)) and two receivers (which receive the signals r \(/) and riU)). In that case, the signals 




received are described by 


ri(f) =Si(t-Z n )+S2(t-Tl2) ( 1 ) 

and 

n{t) = 51 (r- T21) + S2{t- T22), (2) 

where x mn is the delay between transmitter n and receiver m. In the frequency domain, (1) 
and (2) become 

Ri(f) = +S 2 (f)e~R ° Tl2 (3) 

and 

Ri{f) = 5i(/)e-^ T21 +5 2 (/)e”^ T22 , (4) 

where S n (f ) and R m (f) are the frequency-domain representations of ,s„(/) and r m (t), re¬ 
spectively, and <u = 27 t/. Generalizing to a system of m receivers and n transmitters gives 


e ~j(OT 11 . 

. . g j^^ln 


' Si (/) ' 


' *1 (/)" 


. . g j^^mn 


. s " l f '>. 


_ Rm(f) _ 


With an understanding of the delays x mn between each transmitter and receiver, it is possible 
to solve (5) for the transmitted signals S\(f).. .S n (f). It is important to recognize that 
each element of (5) is a function of frequency and must be solved for each frequency. 
Application of the fast Fourier transform (FFT) to the sampled time-domain signals yields 
signal descriptions in the discrete frequency domain and facilitates this solution. 


OFDM Signal Recovery 

A basic system of two receivers collecting the signals of two OFDM transmitters is ad¬ 
equate to demonstrate the potential recovery performance. The transmitter and receiver 
geometry is depicted in Figure 1, and the corresponding distances used to compute the sig¬ 
nal propagation delays are provided in Table 1. The OFDM signal parameters for each of 
the transmitters are summarized in Table 2. 

xviii 



Table 1. Transmitter and receiver geometries. (In kilometers.) 


Element x y alt 

txT o o (T 

Tx 2 0 30 0 

Rx 1 146.2 0 400 

Rx 2 146.2 50 400 



Front Top 


Figure 1. Overview of test geometry. 

Table 2. OFDM test signal parameter summary. 

Value 

Subcarriers 6 

Modulation Scheme QPSK 
Bandwidth (kHz) 600 

Symbols 16 


In order to quantitatively assess recovery performance, the correlation coefficient 


X(M) 


s*s 

7m 


( 6 ) 


is used to compare the original signal s with the recovered signal s. The scalar product 
of two signals is a means of assessing the similarity of two signals [4]. The correlation 
coefficient normali z es the scalar product by signal powers, such that %§ varies from -1 
to 1, where a value of 1 indicates proportional signals, 0 indicates orthogonality, and -1 
indicates proportional signals with different sign convention [5]. 


The correlation coefficients of the recovered signals are provided in Table 3 and indicate 
that the recovered signals are very closely matched to the original signals. A graphical 
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comparison of a subset of the time-domain samples of the original and recovered signals 1 
is provided in Figure 2, and a frequency-domain comparison, a subset of the total signal 
spectrum, is provided in Figure 3. Both the time and frequency perspectives indicate excel¬ 
lent recovery performance. The frequency-domain perspective does indicate a discrepancy 
at the f — 0 point, which is expected as the delay matrix in (5) is singular at that point. 

Recovered signal error is computed using 

S 5 =S-S (7) 

in the frequency domain, where S is the vector of original signal samples and S is a vector 
of the recovered signal samples. Similarly, time-domain error is computed using 

s 5 = s-s, (8) 

where s and s are the time-domain samples of the original and recovered signals. 
Frequency- and time-domain recovery error plots are provided in Figure 4. Given the 
known singularity at / = 0, that point is omitted from the frequency-domain error plot. 
Both error plots indicate good recovery performance, with frequency-domain errors on the 
order of computer rounding error. The steady-state error in Figure 4b is likely a result of 
the singularity at / = 0. 

Conclusion 

The above results show good correlation between the original and recovered signals. While 
additional work is required to ensure that the recovered signals sufficiently isolate desired 
signals from co-channel interference, the method appears to demonstrate the potential to 
support that isolation in a passive SAR system. 

Table 3. Recovery performance figures. 

Correlation 

Signal 1 0.9995 - 4.6035 x 1(T lf 7 
Signal 2 0.9993 + 3.0037 x 10~ 16 / 


1 All signals are simulated and assumed to be either voltage or current. Units are unspecified in the figure. 


XX 




30 40 

(isecs 



Figure 2. Time domain of transmitted and recovered signals. 
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Figure 3. Frequency domain of transmitted 


and recovered signals. 
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CHAPTER 1: 

Introduction 


Synthetic aperture radar (SAR) imaging was conceived in 1951 by Carl Wiley, and the util¬ 
ity of the technique was quickly recognized. The early work in the field depended upon the 
ability to finely resolve the Doppler characteristics, or changes in the signal return due to the 
relative motion between radar and target, of the received radar signal [1], Initial approaches 
relied upon optical processing techniques and utilized a variety of spherical, conical and 
cylindrical lenses to create the radar image [1]. Advances in modern computing and dig¬ 
ital signal processing in the 1980s facilitated the more modern wavefront reconstruction 
techniques that commonly see use today [2]. SAR technological developments have facil¬ 
itated a proliferation of SAR imaging systems including the Shuttle Imaging Radar (SIR- 
C/X-SAR) and Canada’s RADARSAT mission, which are both examples of SAR satellites 
developed for Remote Sensing missions studying the earth [1]. Regardless of the ultimate 
use, SAR provides a unique means to monitor the earth that is capable of overcoming some 
considerable limitations of the more commonly known electro-optical systems. As an ex¬ 
ample, the study of tropical regions by remote sensing scientists was particularly hampered 
by cloud cover in those tropical regions as even after a decade of imaging with the electro- 
optical and infrared sensors on the LANDSAT satellites, there were still regions of the earth 
that had not been imaged without cloud cover [1]. SAR provides a means to collect useful 
data for those regions despite the near constant cloud coverage. 

1.1 Overview 

Much like the benefits afforded to the remote sensing community, SAR offers a number of 
advantages over electro-optical imaging systems to the military and intelligence commu¬ 
nities. The ability to image a target area when obscured by clouds is of equal importance 
to the military and intelligence communities. The ability to image a target without a de¬ 
pendence on illumination from the sun allows targets to be monitored at all hours of the 
day. 

Unfortunately, and unlike the imaging provided by electro-optical systems, the transmit¬ 
ter required for traditional monostatic SAR systems reveals the imaging process when that 
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signal is detected in the target area, affording the target the opportunity to employ coun¬ 
termeasures to prevent the imaging of sensitive target area information. In addition to 
sheltering equipment, the use of an active jammer when the SAR transmitter is detected 
might also preclude the successful imaging of the target area [3]. 

By exploiting transmissions that already exist in the target area, a passive SAR system mit¬ 
igates some of the risks associated with operations targeting an adversary aware of the col¬ 
lection. Gutierrez del Arroyo [4] highlights a number of potential advantages of a passive 
SAR system, including both covertness and decreased susceptibility to electronic attack. 
The passive exploitation of existing signals necessitates a radar designed to operate with 
its receiver well separated from the transmitter. Radar systems that exhibit this separation 
between the radar transmitter and receiver are known as bistatic radars [5]. 

A variety of signal sources might be considered for use with such a passive SAR sys¬ 
tem. Two commercially developed systems, Lockheed Martin’s Silent Sentry and Thales’ 
Homeland Alerter 100, utilize analog FM-radio transmissions for the detection and track¬ 
ing of aircraft. These are the first two known commercial systems to utilize passive radar 
techniques, and the random nature of the FM transmissions that they utilize as illuminators 
provides a signal processing challenge [4]. Orthogonal frequency-division multiplexing 
(OFDM) is a digital waveform that enjoys wide use in communications applications. It has 
also been proposed and extensively studied as a potential radar waveform, with a number 
of characteristics that make it well suited to the radar application [4], [6]. One of the more 
prevalent uses of the OFDM waveform is in fourth generation cellular communication net¬ 
works including the IEEE 802.16 Wireless Metropolitan Area Networks compliant WiMax 
technology and long term evolution (LTE) [7], [8]. The proliferation of this signal coupled 
with the low range and doppler sidelobes of its ambiguity function makes it an excellent 
source for passive SAR imaging [9], [10]. 

That is not to say, however, that using OFDM-based cellular transmissions is without chal¬ 
lenges. Fundamental to cellular technology is the concept of spatial frequency reuse. In 
order to support a large number of users with a limited amount of available spectrum, the 
cellular systems must make use of the same frequencies multiple times in a given geo¬ 
graphic area. As explained in [11], a given coverage area is subdivided into cells (often 
modeled as hexagons, but much more irregular in reality). See Figure 1.1 for an illustration 


2 



of potential cell layouts. Each cell is assigned a subset of the available frequencies, such 
that no adjacent cells utili z e the same frequencies. A frequency reuse distance is estab¬ 
lished such that co-channel interference does not impact the communication channel, and 
only outside of that distance from a cell may the frequencies be used again. This cellular 
allocation of frequencies facilitates an increased density of transmitters, which in turn al¬ 
lows the cell phones and transmitters to operate at lower power levels. This, in combination 
with line-of-sight limitations, mitigates interference between cells [12]. 



Figure 1.1. Frequency reuse in adjacent cells, after [11]. 

While spatial frequency reuse is a convenient way to maximize the utilization of the avail¬ 
able spectrum, it poses a challenge to an asset attempting to use such a signal as an illu¬ 
minator for a SAR image. The SAR image reconstruction process utilizes a matched filter 
and, therefore, requires a detailed understanding of the signal used as the illuminator for 
the image. It is likely that the signals from surrounding transmitters will also be received, 
presenting a co-channel interference issue that must be mitigated before the SAR image 
can be formed. That mitigation is the focus of this effort. 

1.2 Literature Review and Prior Work 

Levanon introduced the use of multi-carrier phase coded signals (now more commonly 
referred to as OFDM) and presented an analysis of their characteristics as a radar illumi¬ 
nator [10]. Gutierez del Arroyo and Jackson have published a number of papers related to 
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the use of commercial transmitters as illuminators for passive SAR, including the devel¬ 
opment of a bistatic ambiguity function for WiMax emitters [9], collection and processing 
of WiMax transmissions for bistatic imaging [6], and analysis of the impact of WiMax- 
specific signal characteristics [4]. That work assumes that interference from neighboring 
cells is mitigated through network design, which is reasonable at low collection altitudes 
but becomes problematic with collections from higher altitude given the increased line-of- 
sight to numerous cells in the network. 

Blind source separation (BSS) is a field dedicated to the recovery of original signals when 
numerous signals have been combined prior to receipt. There are a number of approaches 
to signal separation that might be applied to this problem. Stone presents an overview 
of one BSS approach know as independent component analysis (ICA) [13]. While some 
work in the BSS field has been with a limited number of receivers [14], the bulk of that 
research is for situations with equal numbers of transmitters and receivers as well as for 
real-valued signals, often from audio sources. Douglas has extended blind separation to 
complex valued signals [15], and many of the leading algorithms for blind source separation 
are collected in [16]. Those techniques often depend upon the non-Gaussianity of the 
original signals, and OFDM signals were not readily separated in initial attempts to apply 
those tools. 

Fason [17] and Holder [18] developed a SAR model in MATLAB that can be used for 
analysis of various parameters that impact SAR performance. The work described in this 
thesis is intended as a precursor to the extension of that model to the bistatic case and 
performance analysis using recovered OFDM-based signals. 

1.3 Outline 

An overview of the separation approach is provided in Chapter 2, including a description of 
the methods used for calculating the necessary signal delays. A number of signal recovery 
performance parameters are proposed, and separation is evaluated with both narrowband 
and wideband signals. The frequency response that results from delayed or shifted ver¬ 
sions of those signals is evaluated. Recovery performance under noisy conditions is also 
evaluated. 

A description of the OFDM signal model used is provided in Chapter 3 along with an 
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analysis of the recovery of a number of OFDM signals with various constituent parameters. 
That work is extended to several collection geometries in Chapter 4, and a means to evaluate 
potential collection geometries prior to signal collection in order to increase the likelihood 
of viable signal recovery is also discussed. 

Results are summarized in Chapter 5, and a number of areas requiring further study are 
identified and discussed. 

An analysis of two different methods of imposing a time delay on a sampled signal is pro¬ 
vided in Appendix A. Figures detailing the results of the signal recovery efforts in Chapters 
3 and 4 are provided in Appendix B. MATLAB functions that implement critical elements 
of the work described throughout this thesis are provided in Appendix C. 
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CHAPTER 2: 
Signal Recovery 


Cellular technologies capitalize on line-of-sight limitations in order to use available spec¬ 
trum efficiently. These line-of-sight limitations are negated when the receiver is elevated, 
with potential to receive multiple signals simultaneously. In many applications, use of those 
signals requires the ability to isolate a single transmission from a multitude of co-channel 
transmitters, often without detailed knowledge of the transmitted signals. 

The main contribution of the work described in this thesis is a method to separate co¬ 
channel signals by leveraging the known positional information of both the transmitters 
and receivers to facilitate the recovery of the desired signal. 


2.1 General Approach 

Given a transmitter sending a signal s(t), a receiver will receive r(t) — s(t — t), where 
T is the delay due to propagation between the transmitter and receiver, a function of the 
distance and propagation rate of the intervening medium, and attenuation and noise effects 
are neglected. Considering a small system of two transmitters and two receivers, the signals 
received are 

r\{t) = sx(t - x n ) + si{t - Tn) ( 2 . 1 ) 

and 

riif) = Si(t-T 2 l)+S 2 (t-T22), (2.2) 

where x mn represents the delay between receiver m and transmitter n. Any differences in 
signal losses across receivers are assumed to be negligible. 

Applying the Fourier Transform to (2.1) and (2.2) gives 

Ri(f) = St(/>“ ; '“ Tl1 +S 2 (f)e-^ 2 (2.3) 

and 

Ri{f) = S\{.f)e- j(a ^ +S 2 (f)e -^ 2 , (2.4) 
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where (0 = 2k f and S n (f) and R m (f ) are the Fourier transforms of the transmitted and 
received signals, respectively. In matrix form, and extended for m receivers and n transmit¬ 
ters, (2.3) and (2.4) become 


e ~jO)T n . 

. . g j®^ln 


"Si(/) ' 


" *1 (/)" 

g jCOTml 



_ s„(f) _ 


. «*(/). 


(2.5) 


Given adequate estimates of the delays T mn , (2.5) can be viewed as the canonical Ax = b 
linear algebra problem for all points in /, where A and b are known, and linear algebra 
techniques can be applied to recover x, which holds the frequency-domain description of 
each transmitted signal. 


Working with sampled data of length N, the discrete Fourier transform (DFT), typically im¬ 
plemented via the fast Fourier transform (FFT), yields N frequency bins, with the frequency 
of each bin given by 

(0 = 2 ( 2 . 6 ) 

in radians per second (rad/sec), where i is the bin number and f s is the sampling frequency, 
in Hertz (Hz). This means that there are N separate instances of (2.5) to solve. It is con¬ 
venient to reshape the problem from N unique matrix equations to a single equation that 
incorporates all N frequencies. In that case, the exponential delay terms are arranged in 
m x n blocks on the diagonal of a sparse matrix and (2.5) becomes 


’ A(l) 

0 


" S(l) ' 


" R(l) " 

0 

A (AO _ 


_ §(A0 _ 




(2.7) 


where 


M fs'Z 1 1 


y27T pr fs^ln 


A(/)= ; ; (2.8) 

g~j- n ~W fs^rnn 

and S and R are column vectors of S(i) = [Si(i)... ^(i)] 7 and R(i) = [Ri(i) ■ ■ .R m (i)] T . 
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For convenient reference, (2.7) is summarized as 


AS = R (2.9) 

where A, written without an index, refers to the complete set of A(/) blocks arranged on 
the diagonal, as in (2.7), while S and R are the complete columns of S and R, respectively. 


2.2 Geometry and Delay Overview 


L ry (2)~ 

Lty( 2 )- 




* 





3 3 



(a) Top view 




Figure 2.1. Geometry determines the signal delay. 


The delay z in the above equations is the propagation delay, which is a function of the 
geometric relationship of the transmitters and receivers. An example configuration of two 
transmitters and two receivers is depicted in Figure 2.1, but this is readily extended to 
accommodate additional transmitters and receivers. For simplicity of the distance calcula¬ 
tion, this configuration assumes that the receivers travel on parallel paths, though there is 
no requirement for such a geometry. The configuration measures transmitter and receiver 
location from an origin point and records those locations either in column vectors L rx and 
Lyy for receivers or L tx and L ty for transmitters. The perpendicular slant distance between 
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the receiver path and transmitter d Xmn is computed first, using the law of cosines as 


d Xmn = \j(alt m + R e ) 2 +R 2 e ~ 2(alt m + R e )R e cos (0), 


( 2 . 10 ) 


where 


0 = 


Lrxipi) 


( 2 . 11 ) 


is the interior angle between the receiver path and the transmitter location, as depicted in 
Figure 2.1b. With the perpendicular distance, the angle 7 is defined by 


tan 7 = 


Lry(rri) Liy[l{) 


( 2 . 12 ) 


and the desired distance d mn can be determined using 


cos 7 = 


dr 




(2.13) 


Substituting (2.12) into (2.13), the direct distance between the transmitter and receiver is 
provided by 


dmn — 


dr 


cos 


tan 


— 1 f Lry(m) L;y (// ' 


(2.14) 


The propagation delay is then calculated using z mn = d mn /c, where c is the speed of light. 
A MATLAB implementation of this approach is provided in Appendix C.l. 


A straightforward solution to the canonical Ax = b problem requires an invertible A ma¬ 
trix. In this scenario, that requires both a cooperative geometry and an equal number of 
transmitters and receivers. Geometry can be evaluated in advance, which should enable 
the identification and avoidance of unfavorable collection geometries. The approach is first 
evaluated in the m = n case, but the proliferation of transmitters coupled with the limited 
receiver assets makes the m < n condition worthy of additional consideration. 


2.3 Narrowband Evaluation 

For m = n, the transmitted signals can be recovered as long as the matrix of phase delays 
is invertible. Such a geometry was used to evaluate the approach using narrowband signals 
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in an effort to evaluate the validity of the approach. 


2.3.1 Test Geometry 

The values used in the following analysis are summarized in Table 2.1. The receiver loca¬ 
tions along the .r-axis were calculated to give a look angle a of 20° from nadir. Locations 
along the y-axis were selected to avoid complementary delays. While such delays would 
have resulted in a delay matrix with orthogonal eigenvectors and simplified the analysis 
somewhat, that is not necessarily representative and was, therefore, avoided. 

Table 2.1. Summary of geometry parameters. 


Description 

Parameter 

Value 

Number of transmitters 

n 

2 

Number of receivers 

m 

2 

Receiver altitude 

alt 

400 km 

Earth radius 

Re 

6371 km 

Look angle from nadir 

a 

20 ° 

Transmitter location (x) 

Ltx 

[0 0]km 

Transmitter location (y) 

Lfy 

[ 0 30 ] km 

Receiver location (x) 

Lrx 

[ 146.2 146.2] km 

Receiver location (y) 

Lry 

[ 0 50] km 


The geometry calculations, as described in Section 2.2, and construction of the A matrix 
defined in (2.8) were implemented in a MATLAB functions delays_f rom_laydown() and 
block_deltas_f rom_delays (), which are provided in Appendix C.l and Appendix C.2. 

2.3.2 Signal Generation and Recovery 

Once the delays between transmitter and receiver pairs T mn are established, it is possi¬ 
ble to generate the test signals as both transmitted and received. The MATLAB function 
calc_signals_nb() is used to generate both the transmitted and received signals and is 
included in Appendix C.3. For the narrowband tests, signals were generated of the form 

s n {t) =A M cos(a) M r) (2.15) 

where A n is the signal amplitude for the nth transmitter (A„ = 1 for all of the cases discussed 
below) and oo n = 2k f„ is the corresponding signal frequency. For the first test, signals were 
selected that consisted solely of a single fundamental frequency of the Fourier transform 
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domain - that is, frequencies that conform to integer values of i in (2.6). In that case, 
an integer number of waveforms appears in the time-domain signal, and the frequency- 
domain representation has two non-zero values, corresponding to ±(o n . That restriction 
was not applied to Q)\ in the second narrowband test described below. Specific parameters 
for each test are provided in Table 2.2. The transmitted signals were then sampled at F s , 
and a corresponding frequency-domain set of samples was created using the FFT. 

The sampled received signals can be constructed as the sum of the transmitted signals, 
appropriately delayed, while neglecting any attenuation, and ignoring any noise effects 
temporarily. As the transmitted signals are defined in continuous time, the delayed versions 
can be computed directly for each of the sample points. Each received signal is then 

n 

Tmi)t (2-16) 

i= 1 

and the corresponding frequency-domain signal is computed with the FFT. 


Table 2.2. Signal and sampling parameters, narrowband tests. 



Description 

Parameter 

Value 

Common 

Sampling frequency 

F s 

1.024 MHz 


Number of samples 

N 

32 


Signal magnitude 

An 

1 

First NB Test 

Frequencies (Hz) 

fl 

32,000 



fl 

187,500 

Second NB Test 

Frequencies (Hz) 

h 

47,685 



h 

187,500 


Having generated received signals that incorporate the properly delayed transmissions, it is 
possible to attempt recovery of the transmitted signals by solving (2.7) for S. The MATLAB 
left-division operator is used to determine that solution in the recover_from_delays () 
function, included in Appendix C.4. Solving (2.7) yields the frequency-domain description 
of the transmitted signals. The corresponding time-domain signals are then computed using 
the inverse fast Fourier transform (IFFT). 
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2.3.3 Performance Evaluation 

In addition to a qualitative comparison of the original and recovered signals in both the time 
and frequency domains, some quantitative assessment of the recovered signals is useful. 
Those methods are described here. 


Power Angle 

As the test signals used are real-valued, the recovered signals should share that property. By 
comparing the signal power contained in the real versus the imaginary part of the recovered 
signal, we can compute a power angle, the orientation of the signal power in the complex 
plane. 


For a sampled signal, the power is given by 



(2.17) 


where N is the number of samples in the signal s, which is a column vector, and s* is 
the conjugate transpose of the column, such that s*s computes the inner product of the 
signal [19]. Signal power is composed of real and imaginary power. The angle of the 
power vector in the complex plane is given by 


Q p = tan 


( 3{s} r 3{s} \ 

UismsiJ 


(2.18) 


where the real and imaginary parts of the signal are denoted with the 91 {} and 3{} opera¬ 
tors, respectively, and the length factors (1/A) from (2.17) cancel. A Q p or power angle (as 
referenced below) of zero indicates that the signal power is completely real. Non-zero val¬ 
ues indicate that the recovery process resulted in a signal that has some amount of reactive 
power, which is unexpected, and serves as a first check of the recovered signal. 


Error Figure 

An understanding of the differences between the original and recovered signal is also a 
useful indicator of performance. Such an understanding is given by the energy or power of 
the difference between the two signals, particularly when normalized by the energy of the 
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original signal. This normalized error energy or error figure is given by 

Error Figure = ( ^ S),( ^ S) (2.19) 

s*s 

where s represents the recovered signal and s the original signal [19]. Zero error figure 
indicates that there is no energy in (s — s) — that is, there is no energy in the recovered 
signal that was not in the original signal, which is the ideal case. 

Scalar Product or Correlation Coefficient 

The similarity of two signals can be assessed by computing the scalar product, which is 
maximal when the signals are proportional [12]. Schwarz’s inequality relates the scalar 
product to the signal powers, such that, in the case of transmitted and recovered signals, s 
and s, 

s*s < (2.20) 

where the power of the discrete time signal is given by (2.17). To facilitate comparison of 
signals with different power levels, it is desirable to normalize to unit power signals, such 
that 

*(§,s) = -=. (2.21) 

The correlation coefficient can vary from -1 to 1, where a value of 1 indicates proportional 
signals, 0 indicates the signals are orthogonal, and -1 indicates that the signals are propor¬ 
tional but with differing sign convention [20]. 

Carlson [12] identifies that the cross-correlation function is the scalar product of two sig¬ 
nals with the second signal delayed by T, and the cross-correlation function measures the 
similarity of the two signals as a function of that delay. For purposes of this effort, we 
are interested in the correlation coefficient when T = 0, which corresponds to the scalar 
product. Signals are described as well correlated when the normalized scalar product ap¬ 
proaches one, and correlation is assumed to be at r = 0. 

Subspace Angle 

Taking a page from linear algebra, we imagine the sampled signals as vectors in N- 
dimensional space. Each vector describes a subspace contained within that A-dimensional 
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space. In that context, identical signals are parallel to each other, pointing in the same 
direction in that iV-dimensional space. The angle between two such vectors is given by 


6 = cos 1 



( 2 . 22 ) 


where || • || denotes the function that computes the length or 2-norm of a vector. A 0 = 0 
result indicates parallel vectors while a 0 — 90° indicates completely orthogonal vectors. 

Recognizing that P s — ||s|| 2 /A = s*s/N shows that subspace angle is the inverse cosine of 
the normalized scalar product. Throughout this document, the subspace angle is reported 
in degrees as an alternative perspective to the scalar product. 


2.3.4 Results 





- * 1 ■ o 


-21 - 1 - 1 - 1 - 1 - 1 - 

-600 -400 -200 0 200 400 600 

Frequency (kHz) 


(a) First signal (b) Second signal 

Figure 2.2. Frequency-domain plot of narrowband signals as transmitted and recovered, f\ = 32 
kHz. 


Table 2.3. Performance figures for f\ = 32 kHz. 


Power Angle (°) 
Error Figure 
Subspace Angle (°) 
Correlation 


Signal 1 
5.1252 x10 29 

1.5111 x 10~ 26 
0 
1 


Signal 2 

5.3304x10 ' M) 

1.112 x 10- 26 +1.7106 x 10~ 43 z 
1.8149 x 10~ 14 - 1.8149 x 10 14 z 

1 + 1.0033 x 10~ 31 / 


The results of the first narrowband test indicate a very strong agreement between the trans¬ 
mitted signals and the numerically recovered transmitted signals. The frequency-domain 
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plots 2 (including both the magnitude and phase of each) of the transmitted and recovered 
signals for the first test, with f\ = 32 kHz, are provided in Figure 2.2. For clarity, the 
phase plots only include data points for those frequencies whose magnitude is at least 1% 
of the maximum value in the magnitude plot. Both recovered signals demonstrate good 
agreement with the transmitted signal, including very small phase error for the non-zero 
frequency components of the signal. The corresponding performance figures for this test 
are provided in Table 2.3. The power angle results show that each of the recovered sig¬ 
nals is almost entirely real, with imaginary parts on the order of computer rounding error, 
^machine- Similarly, this test yields very small error figures, indicating that there is very lim¬ 
ited power difference between the two signals. Lastly, the angle between these two vectors, 
when imagined in TV-dimensional space, is very nearly zero, indicating parallel vectors, 
which is the ideal case. Overall, this indicates that the technique has some promise, and a 
good correlation coefficient in the time domain is expected. 





T® 

¥?* 




o 

1,4 4A° 0 i® °46 


-5 - 1 ——©- 1 - 1 - 1 - 1 - 1 - 

0 5 10 15 20 25 30 35 

psecs 


(a) First test (b) Second test 

Figure 2.3. Narrowband tests, transmitted and recovered signals. 


The time-domain results for the first narrowband test are provided in Figure 2.3a. As 
expected, the recovered signal is closely matched to the original signal. Similarly, the 
time-domain signal results from the second test, using f\ — 47,685 Hz are provided in 
Figure 2.3b. That case shows considerable disagreement between the original and recov¬ 
ered signals. Examination of the frequency-domain plots, provided in Figure 2.4, provide 
insight into the problem. 

“All signals throughout the document are simulated and assumed to be either voltage or current. Units are 
left unspecified in all figures. 
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(a) First signal (b) Second signal 

Figure 2.4. Frequency-domain plot of narrowband signals as transmitted and recovered, f\ = 
47,685 Hz. 


Table 2.4. Performance figures for f\ = 47,685 Hz. 


Signal 1 

Power Angle (°) 0.00213 

Error Figure 10.34 

Subspace Angle (°) 110.5 — 0.01644/ 

Correlation —0.3497 + 0.0002688/ 


Signal 2 

0.0003375 

6.676 

72.01 +2.3432 x lO^z- 
OJOSS - 3.8898 x 10~ 19 / 


Unlike the first test, in the second test, f\ is transmitted as a sinusoid that maps onto many 
DFT indices, since f\ yf n/T samp (where is n is an integer), and this transmitted signal 
demonstrates a smooth phase throughout most of the spectrum. (Strum refers to this as 
leakage [21].) The recovered signal lacks this smooth phase, instead yielding a much more 
erratic phase. While the primary frequency of the second signal appears in the recovered 
version (though with reduced magnitude), there are a number of other frequency compo¬ 
nents present in the recovered version that are not present in the transmitted signal. The 
performance figures, provided in Table 2.4, also indicate poor recovery performance. While 
the power angle remains low (indicating real valued signals), the error figure, subspace an¬ 
gle, and correlation all indicate little agreement between the original and recovered signals. 


2.4 Frequency Response Analysis 

The approach outlined in Section 2.1 depends upon the fact that a known delay in the time 
domain translates to a predictable phase shift in the frequency domain and that the Fourier 
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Transform allows for transition between those domains. While these are known properties 
of the Fourier Transform, the results described above indicate that there are likely some 
details of those properties that were not adequately considered. 


2.4.1 Frequency Response Definition 


In order to study the necessary relationship, the delay was considered as a frequency re¬ 
sponse, such that 


H{co) 


= -jcor 


(2.23) 


where ^[-\ represents the Fourier transform operation. By computing the FFT of a signal 
and a delayed copy, it is possible to compare the resulting frequency response with the 
expected e~-’ cot result. As the approach outlined in Section 2.1 is focused on the phase 
relationships, the frequency response magnitude is less important. Any deviations from 
unit magnitude are canceled by the IFFT operation. The phase relationship is the critical 
component, as without a linear phase shift over varying frequency, the approach described 
in Section 2.1 breaks down quickly. For all of the tests below, the delay is T = 2 T samp , 
where T samp is the sampling interval. 


2.4.2 Fourier Transform Relationships 

The delay property applied in (2.3) and (2.4) is a property of the continuous Fourier trans¬ 
form (CFT). For sampled data, a similar property exists for the discrete time Fourier trans¬ 
form (DTFT). When the sampled data is A-pcriodic, the DTFT can be simplified to the 
DFT, which the FFT implements efficiently and also shares a shifting property that yields a 
linear phase shift with frequency. The important characteristic of the DFT, however, is that 
the original samples must be /V-periodic in order for it to be accurately applied in place of 
the DTFT. Ultimately, this requirement for A-periodic data means that the shifting property 
of the DFT is a circular shift [21]. For data that is not ./V-periodic, this difference yields a 
different phase response for the frequency response, and without somehow accounting for 
this difference, the above approach breaks down. These differences are illustrated below. 


2.4.3 Narrowband Cases 

A narrowband signal of the form described by (2.15) was generated with 32 sampled points. 
The magnitude and phase of the frequency response are provided in Figure 2.5 along with 


18 



a depiction of the signals with two different delay types. The first delay, labeled as “Con¬ 
tinuous time shift” was generated by sampling a delayed continuous time signal, while the 
second delay, labeled as “Circular shift” was computed with a circular shift of the original 
signal. Both signals were delayed by two sample intervals. Frequency responses for both 
delay types as well as the ideal response (e~ 70JT ) are provided in Figure 2.5a. Note that 
while not required in this case, later cases involve magnitude deviations from one on the 
order of machine error. For consistency, the magnitude plot is plotted as \H(co) \ — 1 to 
clearly depict those small deviations. 



(a) Frequency response (b) Delay comparison 

Figure 2.5. Delay analysis of the f\ = 32 kHz signal. 


Recall from Section 2.3.4 that this signal provided excellent recovery performance. While 
the phase plot does not show good correlation between either of the frequency responses 
and the ideal response, it is important to recognize that in the frequency domain, ideally 
this signal has zero magnitude for all but two frequencies. When the frequency response 
is computed, both the dividend and divisor are approaching £ mac hi ne f° r all but these two 
frequencies, and the result shows considerable variation. The two frequencies with non¬ 
zero magnitude are marked with small squares, and each of these data points show good 
agreement across the two delay types and the expected frequency response. The time- 
domain view of the two delay types shows they yield similar results, as at this frequency 
and sample rate, the signal is fV-periodic. 

The second narrowband case, as shown previously, has a markedly poorer recovery perfor- 
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(a) Frequency response (b) Delay comparison 

Figure 2.6. Delay analysis of the f\ = 47.7 kHz signal. 


mance, and the frequency response analysis of that case, provided in Figure 2.6, illustrates 
the problem underlying that performance. As the time-domain signal in this case is not 
composed of a single component in the discrete frequency domain, the magnitude of the 
various frequencies is non-zero and, therefore, yields a much less erratic frequency re¬ 
sponse for both types of shifting. While the circular shift of the signal does appear well 
correlated with the ideal case, the signal that was shifted in the continuous time domain 
prior to sampling has a radically different frequency response, particularly with respect 
to the phase angle. The difference results from the fact that the signal is not TV-periodic, 
which is clear from the difference in the shifted signals presented in Figure 2.6b. There, 
two samples (corresponding to the t = 2 T samp delay) of the discrete (or circular) shift are 
not correlated with the corresponding samples of the continuous time shifted signal. The 
delayed signals (as received) were generated using a time shift of a function defined in 
continuous time. A correct solution for S from (2.7) depends upon a phase shift that would 
have resulted from a circular shift vice a delay of the continuous time signal. An accurate 
recovery requires that the effects of the continuous time delay closely resemble those of a 
circular shift. 


2.4.4 Wideband Cases 

For comparison to the narrowband case, a wideband test case was also conducted, using 
an additive white Gaussian noise (AWGN) signal. To facilitate a more direct comparison 
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between this and the narrowband case, this signal was also generated using 32 sampled 
values. As AWGN is not a time-dependent signal, the continuous time delay was created by 
prepending two samples of AWGN to the front of the signal, while the discrete time delay 
was computed using a circular shift, as in the narrowband case above. Using the same 
format as above, we provide transfer functions for both delay types as well as the ideal 
response (e~ JbJr ) in Figure 2.7a. The time-domain plot of the three signals is provided in 
Figure 2.7b, and as in the second narrowband case, also exhibits two points in the circularly 
shifted signal that do not match those in the continuous time delayed signal. Despite that 
difference, the phase plot of Figure 2.7a shows that while the continuous time shift line does 
exhibit a few points of departure from the ideal and circular shift cases, there is generally 
much better agreement between the continuous time shift curve and the ideal case than in 
the narrowband cases, with the overall trend of the curve following that of the ideal curve. 
This indicates that a wider bandwidth signal might adequately support the signal recovery 
approach described above. 



(a) Frequency response (b) Delay comparison 

Figure 2.7. Delay analysis of the AWGN signal, N = 32. 


The same test was repeated with additional samples to determine if the performance was 
a function of the number of samples of the AWGN signal. With additional samples, the 
number of samples that differ between the continuous delay and circular shift is a smaller 
fraction of the overall signal. An AWGN signal was generated using N = 2048, and the 
plots of the frequency response performance are provided in Figure 2.8a while Figure 2.8b 
depicts the first 32 samples of the complete time-domain signal in order to more clearly 
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show the two different points and the overall two-sample delay of the signals. While there is 
some phase noise present in the continuous time-shifted signal, it is important to recognize 
that this signal is plotted over the [—tt, H-7t] interval which, in some cases, causes what 
might otherwise be a smaller error to appear more significant. Aside from this artifact, the 
continuous time plot does appear more closely aligned with the ideal case than any of the 
previous frequency response tests. This improvement is likely due to the larger number of 
samples. With more samples, the two points that are impacted by the shift differences are a 
much smaller fraction of the total signal and have little impact on the frequency response. 



(a) Frequency response (b) Delay comparison 

Figure 2.8. Delay analysis of the AWGN signal, N = 2048. 


2.4.5 Zero Padding 

To further improve the similarity between the delay that occurs in continuous time and the 
circular shift that is required to realize a e~^ co ' c phase shift in the frequency domain of the 
sampled signal, zero padding and windowing can be applied to the collected signal. 

Intuitively, if a signal is transmitted for a limited period of time, a receiver receives no signal 
while the leading edge of the signal is propagating, and the signal arrives only after the 
propagation delay T. The causal nature of the linear time invariant (LTI) channel requires 
that the samples of the signal arrive in the same order in which they were transmitted. 
Without zero padding, the circular shift associated with that delay in the discrete time 
domain results in some fraction of the samples from the trailing edge of the signal shifting 
to the leading edge, which is not representative of the actual signal received after the delay. 
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By padding a time-limited signal with trailing zeros, the propagation delay results in some 
fraction of those zeros moving to the start of the signal rather than some portion of the 
trailing edge of the original signal, preserving the causal nature of the transmission and 
propagation delay. 

Mathematically, for a sampled signal 

x[n] —x{nT ) forn e [0, AA — 1] (2.24) 


of duration N and integer sample numbers n. zero padding yields a signal 


x[n] =jc(«r)/[ 0 ^_i](n) for n e [0,iV — 1], 


(2.25) 


where 


h[n] 


1, if n G A 
0 , otherwise 


(2.26) 


is known as the indicator function and N is the new length of the signal. Adequate zero 
padding requires 

(2.27) 

-*samp 

where T is the duration of the signal, x,j is the propagation delay and r samp is the sampling 
interval. 


The proceeding evaluation of shifting types was repeated with zero padded AWGN signals 
for signal lengths of N — 32 and N = 2048. The signal delay evaluated in this case is 
ideal, such that it is an integer number of samples, the transmitted signal is time-limited, 
and the delayed signal incorporates the complete transmitted signal. The comparison of 
the continuous time shift and circular shift are provided in Figure 2.9 for the N — 32 case 
and Figure 2.10 for the N = 2048 case. Regardless of signal length, the circular shift and 
continuous time delay are nearly identical, which demonstrates that under ideal conditions, 
the continuous time delay that results from signal propagation can be removed through 
either a circular shift in the discrete time domain or by multiplication of a phase shift in the 
frequency domain. 
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(a) Frequency response (b) Delay comparison 

Figure 2.9. Delay analysis of the AWGN signal, N = 32, with zero-padding. 



(a) Frequency response (b) Delay comparison 

Figure 2.10. Delay analysis of the AWGN signal, N = 2048, with zero-padding. 


2.4.6 Phase Error Calculation 


Rather than rely simply upon a qualitative assessment, it is preferable to evaluate the phase 
errors described above quantitatively. For both the continuously shifted and circularly 
shifted signals, the normalized phase error is given by 


Phase Error = 


ZH(co) — ZH(co)\\ 


( 2 . 28 ) 
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where ZH(co) is the phase angle of the frequency response for the shifted signal (either 
continuous time or circularly shifted) and ZH(co) is the phase angle of the ideal transfer 
function, Ze~ j0IT = — ft)r. Ideally, this phase error is zero. The phase angle results are 
presented in Table 2.5. Recall that in the narrowband case with f — 32 kHz, both the 
delay approaches had considerable error as a result of the numerous frequencies with near 
zero response. Other than that case, the error for the circular shift approach is roughly 
£machine , which confirms the graphical assessment above. For the continuous time delayed 
signals, in the narrowband cases, the phase error is approximately one, which indicates 
that the phase error is of the same magnitude as the ideal phase response itself. Some 
improvement is provided in the wideband, N — 32 case, while the wideband, N = 2048 case 
shows a considerable improvement, validating the use of additional points with an AWGN 
signal to mitigate the differences between the continuous and discrete time shifts. Once the 
signals are zero padded, the continuous time delay precisely mirrors the performance of 
the circular shift. It is important to recognize that the signals used in the zero padded cases 
were time limited at the transmitter, such that the initial samples of the received signals are 
zeros corresponding to the propagation delay and not samples of the signal that had been 
transmitted before t — 0 at the transmitter. 


Table 2.5. Frequency response test normalized phase errors. 

Circular Continuous 


Narrowband, f — 32 kHz 
Narrowband, / = 47.7 kHz 
Wideband, N — 32 
Wideband, N = 2048 
Wideband, N = 32, zero-padded 
Wideband, N — 2048, zero-padded 


0.04727 

1.054 

1.5254 x 10' 16 

1.017 

3.4068 x 10~ 17 

0.1097 

7.6214 x 10' 17 

0.007304 

1.7453 x 10' 16 

1.7453 x 10 

1.7632 x 10' 16 

1.7632 x 10 


2.5 Wideband Test 

Having identified a signal that possesses the necessary phase relationship between the orig¬ 
inal signal and one that has been shifted in continuous time, we can evaluate the ability to 
recover that signal. For continuity, the test geometry remains as described in Section 2.3.1. 
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2.5.1 Signal Generation and Recovery 


Sampled baseband AWGN signals were generated for each transmitter. The signals arrive 
at the different receivers with delays dependent upon the transmitter and geometry. As 
those delays are based on the specific geometry of the transmitters and receivers, they are 
not an integer multiple of the sampling interval, and a straightforward circular shift of the 
discrete time signals is inadequate to incorporate this precise delay. Two methods were 
considered to incorporate the appropriate delay into the transmitted signals — delay via 
phase shift in the frequency domain and interpolation of the sampled time-domain AWGN 
signals. A comparison of those methods is provided in Appendix A, and based upon that 
evaluation, the delay was applied via phase shift in the frequency domain of the transmitted 
signal given the small error and low computational complexity. 

As discussed in Section 2.4.5, it was necessary to incorporate zero padding in the signals 
prior to computation of the FFT. The minimum padded signal length is given by (2.27), 
with z d = max ( z mn ). In the interest of FFT computational speed, the signal was padded to 
the next largest power of two, such that 

A = 2l to82 ( 75 “^)l, (2.29) 

where |~•] denotes the ceiling, or round up to the nearest integer, operation. 

Given the unique copies of the signals with the appropriate delays for each transmitter/re¬ 
ceiver pair, the complete received signal is constructed as 

r m [ n \ = £ i ^- 1 {s[k}e - jk{27t/N){FsTmi) } (2.30) 

i= 1 

where £ is the number of transmitters and n is the time-domain sample index number. 

The MATLAB function calc_signals_wb() used for signal generation is provided in 
Appendix C.5, while the basic signal parameters used are provided in Table 2.6. It is 
important to note that the received signals constructed by calc_signals_wb are delayed 
copies of a time-limited transmission, as opposed to signals of infinite duration that are 
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sampled for a specific period in the receiver. 3 

Table 2.6. Signal and sampling parameters, wideband tests. 

Description Parameter Value 

Sampling frequency F s 1 MHz 

Number of signal samples N 2048 

Padded signal length N 4096 

Desired power spectral density S„ 1 W / Hz 

Recovery of the transmitted signals from those received is again accomplished by solving 
(2.7) for S. 

2.5.2 Results 

Given the large number of points in the transmitted signals, plots similar to those used in the 
narrowband test are less informative when assessing signal recovery performance. Instead, 
subsets of the time and frequency domains are provided in Figure 2.11 and Figure 2.12. 
These show excellent correlation between the two signals, with the exception of the f — 0 
points in both Figure 2.12a and Figure 2.12b. Those points are not recovered when solving 
(2.7) as the equation is singular at A(l), which corresponds to 0 Hz. 4 (See (2.8).) Error 
plots, depicting the difference between the original and recovered signals are presented in 
Figure 2.13. The recovery error in the frequency domain is 

S 5 =S-S, (2.31) 

where S is the vector of frequency-domain samples of the original signal and S is a vector of 
the frequency-domain samples obtained by the recovery process. Note that in the frequency 
domain, the error at / = 0 has been omitted from the plot due to the known discrepancy of 
the recovery process. The recovery error in the time domain is 

= s — s, (2.32) 

where s and s are the vectors of the time-domain samples of the original and recovered 
signals, respectively. 

3 See Section 5.2 for further discussion. 

4 In the case of both WiMax and LTE, this point corresponds to the DC subcarrier, which both standards 
specify should be null. The recovery process simply sets this point to zero [7], [22], 
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The performance metrics identified in Section 2.3.3 are provided in Table 2.7. Note the 
small Power Angle, Error Figure and Subspace Angle figures, and the correlation coeffi¬ 
cient very near one. These all indicate good recovery performance. 


Table 2.7. Recovery performance figures for AWGN, N = 2048. 


Power Angle (°) 
Error Figure 
Subspace Angle (°) 
Correlation 


Signal 1 Signal 2 

2.0427 x 1(T 26 2.073 x 1(T 26 

6.0921 x 10“ 5 0.0001202-1.0917 x 10~ 44 / 

0.4472 — 4.263 x 10 _12 z 0.6283 - 1.6754 x l(T 12 z 

1 + 5.8073 x 10 16 / 0.9999 + 3.2064 x 10 16 z 
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Figure 2.11. Wideband test signal, time domain. 


2.6 Performance in the Presence of Noise 

To be a useful technique, the process described in Section 2.1 must demonstrate reasonable 
performance in the presence of noise. To evaluate the performance with noise, the above 
wideband signal test using AWGN signals was repeated with varying levels of AWGN 
noise at the receivers. No filtering was applied to either the AWGN signal or the AWGN 
noise. 

As above, received signals were constructed using (2.30) with x mn determined by the same 
transmitter and receiver geometry described in Section 2.3.1. Noise signals r\ m were con¬ 
structed for each receiver as real valued random samples of a Gaussian distribution and 
scaled by the standard deviation of the samples in order to ensure unit variance of the 
sampled noise signals. 
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(a) First signal (b) Second signal 

Figure 2.12. Frequency domain of wideband test transmitted and recovered signals. 



(a) Frequency domain (b) Time domain 

Figure 2.13. Error of recovered signals. 


At each signal-to-noise ratio (SNR) test point, the necessary noise power level was com¬ 
puted as 


Pr 

1QSNR/10’ 


(2.33) 


where P, is the power in the received signal and the desired SNR level is in decibels (dB). 
The received signal, at the desired SNR, is 


r m — fm + crtj/j 


(2.34) 
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The signal recovery process described above was performed and the correlation coefficient 
between the original and recovered signals ^(s,s) plotted vs SNR in Figure 2.14. Some 
performance degradation is expected, as the recovery process solves A S — R for the signals 
S, while with noise, AS + 77 = R is a more accurate model of the received signals. The 
approach described in Section 2.1 makes no provision for identifying or eliminating noise, 
and, therefore, even under ideal performance conditions, the recovered signals are noisy 
versions of the original signals. To facilitate an understanding of how the recovery process 
compares to what might be considered the ideal case, Figure 2.14 includes curves for the 
correlation coefficient of the original signal with a noisy version of that signal ^(s,s), 
annotated in the legend as “Original”, where 

s m — s m T OY\ m . (2.35) 

As this is intended only to show the impact of noise on the correlation coefficient, there 
was no attempt to include this noise in the propagation model. 

This data indicates that for comparable signal correlation coefficients, the recovery process 
requires roughly 35 dB higher SNR than is necessary if the recovery process is not uti¬ 
lized. This is a considerable performance degradation, and without mitigation, may limit 
the utility of this approach. See Section 5.2 for a discussion of means to both assess this 
performance impact and potential means to mitigate this degradation. Prior to expending 
the effort to improve performance in the presence of noise, the recovery method should be 
tested with the OFDM signals used by fourth generation communication systems. 



Figure 2.14. Recovery performance versus noise level. 
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CHAPTER 3: 

Recovery Using Orthogonal Frequency-Division 
Multiplexing Signals 


In order to be useful as a co-channel interference mitigation approach, it is necessary to 
ensure that the process outlined above provides adequate performance with representative 
signals. OFDM is central to both LTE and WiMax technologies which offer promise as 
illuminators for a passive SAR imaging process. A brief introduction to the OFDM signal 
and an assessment of recovery performance with OFDM signals is presented below. 


3.1 OFDM Signal Generation 

Fundamentally, OFDM can be considered an extension of frequency-division multiplex¬ 
ing (FDM). While FDM suffers from the potential for interference between carriers on 
difference frequencies (due to non-ideal filters), OFDM avoids this with the use of orthog¬ 
onal frequencies, such that for a symbol duration T and subcarrier frequencies j\ [12] 

I cos(2 nfkt) cos(2 nfjt)dt = 0 V j.k where j y k. (3.1) 

Given the trigonometric identity 

cos 0cos0 = ^ [cos(0 — (j>) +cos(0 + 0)], (3.2) 

(3.1) can be rewritten as 

I cos[2m(f k -fj)]+co&[2m(fk + fj)]dt = 0 \/j,k where j y k. (3.3) 

By selecting \ —/& = A/ = 1 /T. the \\ — /) and fk + fj terms are an integer multiple of 
A/, the integrals are zero over the [0, T] interval and the frequencies are orthogonal. 

The signal that represents a single OFDM symbol is assembled from a series of K subcar¬ 
riers. Each subcarrier is a sequence of symbols to be transmitted, either from independent 
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sources or a selection of K symbols to be transmitted simultaneously from a single stream 
of data. Each subcarrier is independently modulated (e.g., by using quadrature phase-shift 
keying (QPSK)) to produce K streams of complex data X k [n ]. For a given OFDM symbol, 
a data element from each subcarrier is further modulated at an orthogonal frequency f k and 
combined with the other modulated subcarrier signals, as in 

K -1 

w[n] = £ X k [n]e j2nfkn n = 0,1,... K - 1. (3.4) 

k =0 

The inverse discrete Fourier transform (IDFT) provides a convenient and elegant means to 
implement this modulation scheme (see [12] for a detailed discussion), such that 

w[n]=lD¥T{X k [n]}. (3.5) 

The total OFDM signal bandwidth is 


B = KAf. 


(3.6) 


3.2 Recovery Performance with OFDM 

Baseband timelimited OFDM signals were generated using code adapted from the work by 
Holder [18]. The parameters for each set of test signals are provided in Table 3.1. (Note 
that the bandwidth listed there is a two-sided bandwidth.) Those signals were delayed 
and combined as discussed in Section 2.5.1, using the transmitter and receiver geometry 
detailed in Table 2.1. For consistency with the preceding discussions, the OFDM signals 
were interpolated and resampled at 1.024 MHz. Relatively low bandwidth OFDM signals 
(as compared to FTE and WiMax) were used for testing to keep the number of samples 
required for the signal delay reasonable at this sampling rate. See Section 5.2 for alternative 
approaches. 

Recovery of the transmitted signals was performed by solving (2.7) for S. Results of the 
recovery process are provided in Table 3.2 through Table 3.5. Note that the OFDM signals 
are not generally real-valued, and a non-zero power angle is not indicative of a recovery 
shortfall. A power angle near 45° indicates comparable real and imaginary parts of the 
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Table 3.1. OFDM test signal parameter summary. 



Set 1 

Set 2 

Set 3 

Set 4 

Subcarriers 

6 

128 

128 

128 

Modulation Scheme 

QPSK 

QPSK 

QPSK 

BPSK 

Bandwidth (kHz) 

600 

600 

1000 

1000 

OFDM Symbols 

16 

4 

4 

10 


signal, which is expected with the OFDM signals. Subsets of the signals from the first 
parameter set depicted in the time and frequency domains are provided in Figure 3.1 and 
Figure 3.2. Qualitatively, these indicate good recovery performance with the recovered sig¬ 
nal plotted at the same points as the original signal. The singularity at / = 0 Hz is apparent 
in the frequency-domain plots. Time- and frequency-domain plots of recovery error from 
the first set of signals are provided in Figure 3.3. The frequency-domain error indicates 
recovery errors on the order of machine error, with a few exceptional frequencies that have 
slightly more error than most. The steady-state error in the time domain is likely the result 
of the singularity at / = 0 Hz. Each set of OFDM parameters demonstrated comparable 
performance, with some minor differences that are likely the result of the random signal 
content rather than a function of one of the OFDM parameters. The complete set of OFDM 
recovery performance plots are provided in Section B.l. This consistent and accurate re¬ 
covery performance suggests that the method is likely functional across a large collection 
of OFDM-based signals, including LTE and WiMax transmissions. It is important to rec¬ 
ognize that this performance was examined in a single arrangement of transmitters and 
receivers. As it is useful to understand the potential impact of collection geometry on 
recovery performance, this is explored further in the following chapter. 


Table 3.2. Recovery performance figures for OFDM, set 1. 


Power Angle (°) 
Error Figure 
Subspace Angle (°) 
Correlation 


Signal 1 
47.39 

0.001054 + 7.7941 x 10 ~ 29 i 
1.86 + 8.1252 x 10~ 13 i 
0.9995 - 4.6035 x 10~ 16 z 


Signal 2 
35.82 

0.001404-2.8107 x 10~ 30 / 
2.147 -4.5929 x 10- 13 ^ 
0.9993 + 3.0037 x 10~ 16 / 
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Table 3.3. Recovery performance figures for OFDM, set 2. 



Signal 1 

Signal 2 

Power Angle (°) 

45.45 

44.38 

Error Figure 

0.0004209 — 1.6937 x 10~ 28 z 

0.0004143 + 2.4529 x 10~ 28 / 

Subspace Angle (°) 

1.176 + 4.1564 x 10“ 13 i 

1.166 — 1.7555 x 10 -13 / 

Correlation 

0.9998 — 1.4883 x 10~ 16 « 

0.9998 + 6.2367 x 10- 17 / 


Table 3.4. Recovery performance figures 
Signal 1 

for OFDM, set 3. 

Signal 2 

Power Angle (°) 

42.56 

47.23 

Error Figure 

0.0005003 + 2.2331 x 10~ 28 z 

0.0004999+1.6457 x lO^ 28 * 

Subspace Angle (°) 

1.282- 1.9555 x 10“ 12 i 

1.281+3.2425 x 10 ~ u i 

Correlation 

0.9997 + 7.6339 x 10 16 f 

0.9998 - 1.2653 x 10- 17 / 


Table 3.5 

. Recovery performance figures for OFDM, set 4. 


Signal 1 

Signal 2 

Power Angle (°) 

45.54 

42.5 

Error Figure 

0.0001 - 2.2446 x lO" 30 / 

0.0001001 - 1.7213 x 10~ 29 z 

Subspace Angle (°) 

0.573 + 1.1413 x 10 -13 / 

0.5732- 1.3908 x 10~ 12 z 

Correlation 

0.9999 - 1.9919 x 10 17 z 

0.9999 +2.4284 x 10~ 16 i 
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(a) First signal (b) Second signal 

Figure 3.2. Frequency domain of OFDM test transmitted and recovered signals. 
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Figure 3.3. 


(b) Time domain 
Error of recovered OFDM signals. 
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CHAPTER 4: 

Collection Geometry Analysis 


The recovery method in Section 2.1 depends upon an invertible A matrix in AS = R, as 
identified in Section 2.2. In addition to an equal number of transmitters and receivers (m = 
n) such that A is square, the collection geometry plays a significant role in the invertibility 
of this matrix, as each element of the A(i) block is of the form e ~ j27cl N~^ Xmn and depends 
upon the propagation delay x mn . A means to assess a given collection geometry and an 
initial analysis of the impact of several such geometries upon the recovery performance is 
provided in this section. 

4.1 Geometry Assessment 

In linear algebra, the sensitivity of x to small perturbations in b when solving the Ax = b 
problem is described by the condition number of the m x m A matrix, 

k(A) = — (4.1) 

O m 

where Oi is the largest singular value and o m is the smallest [23]. (Both are computed 
via the singular value decomposition of A.) A matrix with a small condition number is 
described as well-conditioned, while a matrix with a large condition number is described 
as ill-conditioned. A singular matrix is often indicated by k( A) = [24]. Condition 

numbers that yield acceptable or unacceptable results are very dependent upon the details 
of the specific problem. See [23] or [24] for a more complete discussion of the condition 
number and singular value decomposition. 

In the case of (2.7), the A matrix is a block diagonal matrix consisting of the A (i) blocks 
associated with each frequency ft), as defined in (2.6). As each of these blocks is indepen¬ 
dent of the other blocks in the A matrix, it is possible to assess the condition numbers of 
the blocks independently, for the given frequency associated with that block, rather than 
simply evaluating the condition of the complete A matrix. This assessment provides in¬ 
sight as to the recovery performance at each frequency used in the recovery calculation and 
is provided in conjunction with the frequency-domain recovery error plot for each of the 
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layouts presented in the following section. 


4.2 Layout Overview and Test Signal 

Four transmitter and receiver geometries are evaluated, in addition to the layout used in the 
preceding sections, which is included as a reference layout (Layout 1). The transmitter lo¬ 
cations remain constant throughout the analysis. Layout 2 extends the separation between 
the receivers to 500 km from the previous 50 km. Layout 3 uses this same extended re¬ 
ceiver baseline and increases the altitude of the second receiver from 400 km to 800 km. 
Layout 4 moves the second receiver such that it appears at roughly five degrees above the 
transmitters’ horizon. Those receivers have an altitude of 200 km. Lastly, Layout 5 eval¬ 
uates the possibility of co-locating receivers on a single platform using a 1-meter receiver 
separation. These layout parameters are provided in Table 4.1, with parameters that vary 
from the reference layout indicated in bold font. An illustration (not to scale) of the five 
layouts is provided in Figure 4.1, which includes a red outline of the reference layout. 

To facilitate comparison across test layouts, a single set of OFDM signal parameters was 
utilized for all tests. Those signal parameters (Set 1 from Section 3.2) are provided in 
Table 4.2. Additionally, in order to minimize any variation in the recovery performance 
and facilitate comparison across the various geometries, the content of each transmitter’s 
signal was maintained as a constant across each of the layouts. 

4.3 Results 

Correlation coefficient values for each recovered signal are provided for the various lay¬ 
outs in Table 4.3. The complete set of metrics discussed in Section 2.3.3 are provided in 
Appendix B.2. Each case has a correlation coefficient at or very near one, and with one 
exception, all have imaginary parts on the order of £ mac hine- The only exception is with 
Layout 5, with the two receivers on the same collection platform, and in that case the imag¬ 
inary part is between two and three orders of magnitude larger than £ m achine but remains on 
the order of 10 13 . 

Frequency-domain recovery error (as described in Section 2.5.2) and condition number 
plots are provided for each layout in Figure 4.2 through Figure 4.6. The condition number 
plots provide both an overview of the full spectrum used for recovery, as well as a view of 
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Table 4.1. Transmitter and receiver geometries, by layout. (In kilometers.) 



Element 

X 

y 

alt 

Layout 1 

Tx 1 

0 

0 

0 

(Reference) 

Tx 2 

0 

30 

0 


Rx 1 

146.2 

0 

400 


Rx 2 

146.2 

50 

400 

Layout 2 

Tx 1 

0 

0 

0 


Tx 2 

0 

30 

0 


Rx 1 

146.2 

0 

400 


Rx 2 

146.2 

500 

400 

Layout 3 

Tx 1 

0 

0 

0 


Tx 2 

0 

30 

0 


Rx 1 

146.2 

0 

400 


Rx 2 

146.2 

500 

800 

Layout 4 

Tx 1 

0 

0 

0 


Tx 2 

0 

30 

0 


Rx 1 

146.2 

0 

400 


Rx 2 

1113 

50 

200 

Layout 5 

Tx 1 

0 

0 

0 


Tx 2 

0 

30 

0 


Rx 1 

146.2 

50.001 

400 


Rx 2 

146.2 

50 

400 


Table 4.2. Geometry analysis test signal parameter summary. 

Value 

Subcarriers 6 

Modulation Scheme QPSK 
Bandwidth (kHz) 600 

Symbols 16 


a subset of the lower frequencies, depicted on a logarithmic y-axis. (The area contained in 
the subset view is indicated in the plot of the full spectrum with a red box.) This provides 
a better understanding of the periodic nature of the condition number that occurs with 
some geometries. As the condition number of the A(0) block is known to be infinite, it 
has been omitted from the plots. In each layout, there is a clear correlation between the 
condition number at a given frequency and the frequency-domain recovery error depicted 
in the corresponding figure. It is also worth noting that the largest condition number (which 
occurs in the Layout 5 case) is also correlated with the largest recovery error. Even in that 
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case, with a condition number on the order of 10 7 , the frequency-domain recovery error 
only approaches 6 x 10 l4 , and the correlation coefficients for that layout remain close to 
one, indicating excellent recovery performance. 

The complete set of signal and error plots for each of the layouts is provided in Ap¬ 
pendix B.2. Note that Layouts 2, 3, and 4 involve larger propagation delays than Layouts 1 
and 5 and require additional zero padding in the time domain. That additional zero padding 
results in samples at a closer interval in the frequency domain, which is apparent in the 
associated figures. 

Table 4.3. Correlation coefficient summary by layout number. 

Signal 1 Signal 2 

Layout 1 0.9987 - 2.6904 x 10 1V / 0.9991 - 2.496 x 10 lfo z 

Layout 2 0.9994+ 1.9458 x 10 16 / 0.9995 - 7.809 x 10 16 / 

Layout 3 0.9994 - 6.9239 x 10 17 / 0.9995 - 5.1021 x lO’ 16 / 

Layout 4 0.9994-2.9356 x 10 17 / 0.9995+ 1.2612 x 10~ 16 / 

Layout 5 0.9987 - 2.2326 x 10 ~ n i 0.9991 + 6.5599 x 10~ 13 / 
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Figure 4.1. Overview of geometry test laydowns. 
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Figure 4.2. Geometry induced recovery errors, layout 1. 
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Figure 4.3. Geometry induced recovery errors, layout 2. 
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Figure 4.4. Geometry induced recovery errors, layout 3. 
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Figure 4.5. Geometry induced recovery errors, layout 4. 
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(a) Frequency-domain error (b) Condition number 

Figure 4.6. Geometry induced recovery errors, layout 5. 


44 














CHAPTER 5: 
Conclusion 


An approach to signal recovery based upon the known propagation delays between a set 
of known transmitters and receivers was described in this thesis. While results with nar¬ 
rowband signals do not offer much initial promise for the approach, by considering the 
frequency response associated with delayed signals and the need for the delay in contin¬ 
uous time to resemble the circular shift associated with discrete time signals, admirable 
performance was demonstrated with the wideband signals typical of fourth generation cel¬ 
lular transmissions. 


5.1 Summary of Results and Identified Limitations 

Examination of recovery results with several OFDM signals as well as in a variety of 
geometries indicates that the proposed method is worth further consideration and study. 
The recovery approach produced correlation coefficients approaching one for the various 
OFDM signal parameters evaluated as well as the variety of collection geometries. The 
use of the condition number to evaluate prospective collection geometries offers insight 
into potential errors in the frequency domain of the recovered signals and the viability of a 
potential collection from a geometry perspective. 

A number of limitations exist that must be addressed before the approach is viable for real- 
world signal recovery. Performance in the presence of noise is a considerable limitation. 
The above results indicate the need for an additional 35 dB of SNR, a considerable margin 
under any realistic circumstances. The known singularity at / = 0 Hz may also be a limi¬ 
tation if applied to a signal that has a non-zero DC subcarrier. Computational complexity 
must also be considered, given the potential size of the A matrix. These limitations and 
potential solutions are further described in the following section. 

5.2 Future Work 

There are a number of areas that merit additional study, either to address identified short¬ 
comings of the approach as it currently exists or to extend the work in order to ensure 
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adequate performance in the intended application. 


5.2.1 Received Signal Model Improvement 

As discussed in Section 2.5.1, the model currently constructs wideband signals as the sum 
of delayed copies of time-limited transmissions. For small differences in delay times, this 
results in a small fraction of each received signal that only contains signal from a single 
transmitter, as, due to differences in the propagation delays between the transmitter and 
receiver pairs, the signals from the other transmitters have either not yet arrived or have 
completely arrived. For longer propagation delays coupled with shorter signals, this can 
potentially result in the signals arriving at the receiver at distinctly different times. 

While the recovery process does not take advantage of this difference in arrival times, it 
is more likely that a base-station transmitter would exhibit a much larger transmit duty 
cycle and a more accurate model would ensure that the each sample of the received signal 
contained a contribution from each of the transmitters. This would require the generation 
of longer duration signals and the application of a windowing or indicator function such 
as (2.26), where the indicator domain A is limited to a period of time when all signals are 
present in the receiver. 

5.2.2 Computational Complexity 

In (2.7), the A matrix is mN x nN. where m and n are the numbers of receivers and trans¬ 
mitters, respectively. The parameter N is the number of discrete frequencies, which is a 
function of the desired signal duration, sampling rate and amount of zero padding required 
for the given signal delay. Increasing the sample rate in order to accommodate higher 
bandwidth signals or the extending the buffer length when zero padding signals with larger 
propagation delays has the potential to dramatically increase the size of the A matrix and 
the computation time associated with signal recovery as a result. Of note, in this work, 
this limitation required a lower altitude for the off-nadir receiver in Layout 4 of the geom¬ 
etry test. It may be possible to use the difference in propagation time across the various 
transmitter/receiver pairs vice the absolute time, which, in some geometries would offer a 
considerable savings in the amount of zero padding required. 

Some improvement may also be available as a result of the fact that the A matrix is par¬ 
ticularly sparse. Use of MATLAB’s sparse matrix tools should both decrease memory 
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requirements and decrease computation time given that the number of non-zero entries in 
the A matrix grows as O(N) while size grows as 0(N 2 ). Similarly, as currently defined, it 
is possible to break the AS = R equation into a set of A(i)S(i) = R (i) equations for each 
discrete frequency i and solve those independently. This would complicate approaches that 
leverage any potential correlation between adjacent frequencies to improve performance 
with noise or in the m < n case. (See below.) 

5.2.3 Accuracy of Delay Knowledge 

The preceding work is entirely dependent upon and assumes completely accurate knowl¬ 
edge of the delays associated with signal propagation between each transmitter and receiver 
pair and uses double precision variables for the associated computations. There is no error 
included in any transmitter or receiver position data, either when constructing the received 
signals or when solving for the signals as originally transmitted. 

Such levels of accuracy are not available when using real-world transmitters and receivers, 
and it is necessary to evaluate performance under varying accuracy conditions. That analy¬ 
sis should likely include both an analytical evaluation of the pertinent recovery equations to 
quantify the impact of a small change in the expected signal delay on the recovered signals 
as well as a Monte Carlo simulation with varying error levels to confirm that analysis. 

The matter is complicated by the interactions between the different delays in the A matrix. 
If each transmitter and receiver has an associated position error, the geometric arrangement 
of those antennas will have a different effect on the delay between each transmitter and 
receiver. Those various errors are then combined in the A matrix and impact the recovered 
signal. It would be useful to understand the impact of a positional error associated with a 
single transmitter or receiver as well as a set of comparable errors across the complete set 
of transmitters and receivers. 

5.2.4 Noise 

As identified in Section 2.6, noise poses a considerable problem for the recovery process as 
it solves the AS = R equation while AS + rj = R is a more accurate model of the received 
signals. One potential approach to mitigating the impact of the noise on signal recovery 
relies upon the fact that the signals of interest are transmitting digital data in accordance 
with published standards, such as [7]. 
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Those standards include a specific symbol library, and forcing the recovery signal to ad¬ 
here to the library may mitigate some of the impacts of noise. After an initial attempt 
at signal recovery, the recovered signals might be processed as they would by a typical 
commercial receiver, including the detection process, yielding channel symbols [25]. By 
re-encoding those recovered channel symbols, the resulting signal is effectively a clean 
version of the recovered signal (precisely what was transmitted, in the ideal case), and any 
noise-induced differences between the recovered symbol and the re-encoded symbol are 
eliminated. While this does not protect against any difference that effectively results in an 
incorrect channel symbol, it may offer some performance improvement over the existing 
recovery process. While it likely involves a considerable increase in complexity, it may be 
possible to extend this decoding/re-encoding process to leverage any error correction cod¬ 
ing present in the signal-of-interest, further mitigating the impact of noise on the recovered 
signal. 

It is also important to recognize that the condition number serves as an indicator of the 
sensitivity of the recovered signal to small perturbations in the received signals. When 
the condition number is larger, the noise in the received signals will have a greater impact 
on the recovered signal than in those cases with smaller condition numbers. With further 
study, it may be possible to develop an empirical understanding of the cases when a larger 
condition number requires some alternative recovery approach or noise mitigation strategy 
prior to or as part of the recovery process. 

5.2.5 Evaluation of Recovered Signals as Bistatic Illuminators 

Ultimately, the signal recovery process is intended to support bistatic SAR imaging, with 
a recovered signal serving as the illuminator for the imaging process. It is important to 
understand the implications to the imaging process when a recovered signal is utilized vice 
the more typical case employing a known pulse shape. Holder [18] developed a model of a 
monostatic SAR system which could be extended to accommodate bistatic signals and then 
used to evaluate the impact of the recovery process. 

There are two fundamental changes to the existing model necessary to support that evalua¬ 
tion. First, the model needs to be extended to accommodate bistatic imaging, including the 
incorporation of bistatic delay calculations for the layout of targets. In the monostatic case, 
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the round-trip fast-time delay is [2] 


2 a Jx 2 + (y-u) 2 
c 


(5.1) 


where x and y describe the target position, u is the radar’s position along the y-axis, and c is 
the speed of light. In the bistatic case, this delay calculation must account for the location 
of the transmitter and is no longer a round trip time. In that case, the delay is given by 


t 


~ \!\x- p ) 2 + (y - q ) 2 + \fx 2 + {y - u) 2 


(5.2) 


where p and q describe the location of the transmitter. These delay calculations are then 
used to construct the signals received during the collection process. 


The SAR model described by Holder uses Spatial Frequency Interpolation for the image 
reconstruction process, which converts the received signals to a recognizable image. This 
algorithm requires two-dimensional, frequency-domain matched filtering [2] 


S((Q,ku)So(co,k u ), (5.3) 

where S(co,k u ) is the frequency-domain representation of the received SAR signal and 
Sq(g),& m ) is the two-dimensional matched filter frequency response. While Soumekh [2] 
describes the derivation of 5 q( 0),A;„) for the monostatic case, it will be necessary to adapt 
that to the bistatic case in order to create the SAR images required to properly evaluate 
signal recovery performance. The recovered signal is used directly in that matched filtering 
process and is referred to as the reference signal in both [2] and [18]. 

The ability to study the impact of using recovered signals for SAR image formation then al¬ 
lows further examination of the impact of noise and positional errors on the overall system. 
While both of those can be worked somewhat independently, that work ultimately impacts 
the quality of the SAR image, and understanding the impacts of both noise and positional 
error on the image quality is critical to the utility of the approach. 
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5.2.6 Adaptation to m < n Scenarios 

Given the proliferation of LTE and WiMax transmitters, coupled with the fiscal constraints 
and limited resource availability, a signal recovery approach that functions when there are 
more transmitters than receivers is particularly desirable. This is a classical example of an 
underdetermined problem, as given the resulting shape of the A matrix, there is not enough 
information available to yield a single solution to the AS = R matrix equation. Common 
approaches to this underdetermined problem such as seeking the solution that contains the 
most zero elements or that minimizes an error function (often referred to as the method of 
Least Squares) do not seem well suited to this particular problem, as there is no expectation 
that the original waveforms adhered to either of those paradigms. 

By oversampling the received signals, it is possible to develop some redundant informa¬ 
tion in those signals. In the frequency domain, this means that the adjacent samples are no 
longer independent. By incorporating that assumption into the construction of the A ma¬ 
trix, such that the solution for a given point in the frequency-domain representation of the 
signal is dependent upon other surrounding values, it may be possible to develop a better 
approximation for the original signal. 
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APPENDIX A: 
Signal Interpolation 


In order to evaluate the signal isolation approach described in Section 2.1, it is necessary to 
generate a collection of signals that contains a set of original signals at various delays. This 
is straightforward in the narrowband case which uses signals defined over a time limited 
period t of the form 

s n {t) =A„cos(gv) (A.l) 

as the cosine is a continuous function and the required delays can be computed as 
s n (t — T mn ). In the wideband case, which simply uses time-limited AWGN as the signal-of- 
interest, one approach to computing the various delays is to interpolate the original signal 
at the necessary delayed sample times. This interpolation is commonly performed using 
Ideal Low Pass Filter Interpolation. As an alternative, an approach utilizing a phase shift 
in the frequency domain was evaluated. Both are described below with a comparison of 
results of the two approaches. 


A.l Ideal Low Pass Filter Interpolation 

Interpolation allows recovery of information from the original continuous time signal that 
is not immediately available after the sampling process. This approach to interpolation 
utili z es a low pass filter to counter the effects of the sampling process. The sampling 
process and its effects in the time and frequency domain are reviewed below prior to an 
explanation of interpolation using an ideal low pass filter. 


A.1.1 Sampling in the Time Domain 

Ideal sampling of a continuous time signal s(t) can be considered as the product of that 
signal with a train of unit impulses 

—|—oo 

p(t)= Y ( A -2) 
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yielding the sampled signal 


x s (t) = x(t)p(t) 

-\-oo 

— £ x(t)8(t — nT) 

n=—°o 

-\-oo 

= £ x(nT)8(t —nT) (A.3) 

n=—oo 

as a train of pulses with area determined by the signal at each sample interval. For con¬ 
venience, the notation x[n] is used to represent a discrete time equivalent of x(nT). This 
process is illustrated in Figure A.la. Once sampled, the continuous time information is no 
longer directly available — interpolation is required to recover signal information at other 
than the sampled instants. 
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(a) Time domain 
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2/v ./; ./| 0 /, f s 2f s 

(b) Frequency domain 


Figure A.l. Sampling a continuous time signal, after [19]. 
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A.1.2 Frequency-domain Perspective 

In the frequency domain, depicted in Figure A. lb, the impulse train used for sampling is 
also an impulse train, and taking f s — l/T, is described by 

-)-oo 

P(f)=fs £ 8 (f - kf s ) . (A.4) 

k=—oo 

The multiplication of s(t ) and p(t) in the time domain corresponds to convolution in the 
frequency domain, such that 

X s (f)=X(f)*P(f) 

-\-oo 

= fs I X(f)*S(f-kf s ) 

k=—o o 
-\-oo 

= fs £ X(f-kf s ), (A.5) 

k=—o o 

which is a train of the original spectrum repeated at multiples of f s in the frequency domain 
[19], [21]. As long as the sampling rate exceeds that required by the Nyquist sampling 
theorem, aliasing is avoided. 


Uf) 



-2fs -fs -fl 0 /; ./; 2f s 


Figure A.2. Ideal low pass filter. 


A.1.3 Interpolation in the Frequency Domain 

In the absence of aliasing, in order to recreate the original continuous time signal, the 
original frequency spectrum can be obtained by applying an ideal low pass filter to isolate 
the spectrum of the original signal, as depicted in Figure A.2 [12]. The ideal low pass filter 
is described by 

Hilpf(Z) = (4) (A ' 6) 
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where 


n< 


x 


I if |^| <1/2 
0 otherwise 

is known as the rectangle function. Applying the ideal low pass filter to (A.5) gives 


(A.7) 


Y(f) = / W/)W) = Kf s X(f), (A.8) 

and selecting K — 1 / f s gives Y (/) — X(f). The inverse Fourier transform can then be used 
to recover the original time-domain signal. 


A.1.4 Interpolation in the Time Domain 



Figure A.3. Signal interpolation as a sum of scaled and shifted sinc() pulses. 

The multiplication of X s (f) by the ideal low pass filter’s frequency response in (A.8) is 
convolution in the time domain, such that y(t) = /zilfp(0 Substituting (A.3) then 

gives 


y(t) =hi LFP (t)*Y,x[n]8(t-nT) 

n 

= J^x[r]/zilfp(0 * 8 (t-nT) 

n 

= Y d x[n\hi L w(t-nT). (A.9) 

n 
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Given the Fourier transform relationship 


*[1 (^j & 2BKsinc(2Bt), (A.10) 

and choosing B — f s /2 and K — 1 /f s , (A.9) becomes 

y(t) = J^x[n] sinc(f s t- n) (A.ll) 

n 

which is a sum of shifted sinc() pulses, each scaled by a corresponding coefficient from 
the x[n] vector [12]. This is illustrated in Figure A.3, which depicts the constituent sincQ 
functions that are added to construct the x(t) signal. Each sinc() pulse is zero at all mul¬ 
tiples of T other than the peak, and each sample point only receives a contribution from a 
the single sinc() that corresponds to that particular x\n\ value. (Note that all but one of the 
sinc() curves are zero below each of the red circles corresponding to a sample in x[n\.) All 
other points in the x(t) signal are composed of contributions from each of the sincQ pulses. 
MATLAB code to perform this interpolation function is provided in Appendix C.6 as the 
function interp_time .m. As MATLAB works with sampled signals, this code returns the 
signal resampled at new time values, which are specified as an argument to the function. 


A. 1.5 Error Analysis 

As a baseline for later sections, it is useful to understand the amount of error induced by 
this interpolation process. As a challenging case, a signal x(t) consisting of AWGN was 
sampled to create x[n], A sampled version of x(l — t) was generated by interpolating s[n] 
at a delayed point in time, where T is an integer multiple of the sampling interval and a 
significant fraction of the overall signal duration to make the delay apparent in the follow¬ 
ing figures. For comparison, a circular shift corresponding to r was applied to the same 
sampled AWGN signal. Those signals are presented in Figure A.4, along with the differ¬ 
ence between those two signals (plotted on a separate y-axis to emphasize the dramatic 
difference in scale.) The total energy in both shifted signals was computed as 

Ei = x\ XjT sam p, (A. 12) 
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where xj is a column vector of the signal-of-interest with sampling interval r sam p, and X|* 
is the conjugate transpose of the column, such that xfxj computes the inner product of the 
signal. This was compared to the energy in the original x[n\. The differences in energy are 
presented in Table A. 1 and show error on the order of £ m achine> which indicates that shifting 
via interpolation is a viable option. 



0 500 1000 1500 2000 2500 3000 3500 4000 4500 

Time sample 

Figure A.4. Signals and error from shifting via interpolation. 


Table A.l. Energy summary for shifting via interpolation. 

Energy Error (J) 

Circular Shift 4.3368 x 1(T 19 
Interpolation 2.4286 x 10~ 17 


A.2 Fourier Phase Shifting for Interpolation 

In the continuous time domain, the Fourier transform has a property that relates a shift in 
the time domain to a phase shift in the frequency domain [19], such that 

x(t)^X(co) (A.13) 

x(t-t d )&X((o)e- j(0td . (A. 14) 

In the discrete time domain, the discrete Fourier transform has a similar property such 
that a circular shift in the discrete time domain results in a phase shift in the frequency 
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domain [21]. That is, 



(A.15) 


(A. 16) 


where n is the index or sample number, © represents the left circular shift operation, m 
is the integer number of samples to shift by, N is the total number of samples, and k is 
the discrete frequency index. With adequate zero padding, this circular shift resembles the 
time delay from the continuous time domain, where the delay is given by t c i — mT samp and 
Tsamp is the sampling interval. This allows for shifting of signals by integer multiples of the 
sampling interval by applying the FFT, multiplying by the appropriate phase shift in the 
frequency domain, and returning to the time domain with the IFFT. While this convenient 
for signals that are already in the frequency domain, it is likely excessive effort for a signal 
in the time domain, as a circular shift operation is accomplished with a straightforward 
manipulation of the index numbers. 

For comparison to other cases, plots of the delayed signals and associated error between 
this method and the circular shift in the time domain are provided in Figure A.5. As this 
method involves multiplication by complex terms, plots of both the real and imaginary parts 
of the signal are provided. Similarly, the energy difference between the phase shifted and 
circularly shifted signals are provided in Table A.2, which also includes entries for both the 
real and imaginary portions of the error. Each is well below £ ma chine- 
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Figure A.5. Signals and error from integer phase shifts in the frequency domain. 
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Table A.2. Energy summary for integer phase shifts in the frequency domain. 

Energy (J) 

Error 8.9355 x l(T 2y 

Real(Error) 8.934 x 1CT 29 
Imag(Error) L51 x 10” 32 

A.2.1 Fractional Delays 

The opportunity to use non-integer or fractional delays for shifting is appealing as it ac¬ 
complishes an interpolation to get the new set of values of x(t) that are required but not 
available directly in x[n ]. The plots of the signals and error that result from such an attempt 
are provided in Figure A.6. As the circular shift does not accommodate fractional shift 
amounts, an interpolation was used as the baseline signal for comparison and when com¬ 
puting both the error plot and energy errors, which are presented in Table A.3. The error 
levels in this case are much more significant than those with integer shift amounts, and, 
more significantly, the error in the imaginary portion of the signal is an order of magnitude 
larger than that of the real portion. Intuitively, as the process begins with a real signal and 
the intent is to merely delay the signal in time, the result should not have an imaginary 
portion. 


Real part of time domain signals 



Time sample 

Figure A.6. Signals and error from fractional phase shifts in the frequency domain. 

Table A.3. Energy summary for fractional phase shifts in the frequency domain. 

Energy (J) 

Error 6.1483 x KT 7 

Real(Error) 7.2528 x 1(T 8 
Imag(Error) 5.423 x 10 -7 
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Error Source 

The Fourier transform, in both the continuous and discrete forms, has a number of prop¬ 
erties that are helpful when considering this error. In particular, when x(t) is real, the 
frequency-domain description of the same signal X(f) will be conjugate symmetric, such 
that X(—f) — X*(f) [19]. The inverse is also true, such that a signal must be conjugate 
symmetric in the frequency domain in order to be real in the time domain. 

Recall that in the case of discrete signals, the frequency domain is periodic. In that case, 
conjugate symmetry requires that X[n] = X*[N — n] where the signal is periodic over A 
points, n < N, and the first point is n = 0. (In MATLAB, indices begin at n — 1, so in that 
case conjugate symmetry requires X[n] — X*[N + 2 — n\. The remainder of this discussion 
will assume the first index is n = 0.) Two of these points are significant — n — 0 and 
n = A/2. When n = 0, the above relationship requires conjugate symmetry with the point 
at n=N, which exceeds the length of the sampled signal but corresponds to itself when the 
periodic nature of the DFT frequency domain is considered. The same is true for n — N/2, 
which must also be conjugate symmetric with itself, or, more simply, real. 

When shifting a signal by a fraction of the sampling interval, the e ~ jC,n term must be 
conjugate symmetric in order for the resulting signal to be real-valued in the time domain. 
That requires that both the n = 0 and n — N/2 values in that term to be real valued, which 
is not the case when shifting by fractions of a sampling interval. In that case the n = A/2 
value will be complex, and the resulting signal will have an imaginary part in the time 
domain. 


Forcing Conjugate Symmetry 

As a real valued result is expected, the n = A/2 point can be artificially forced to a real 
value by eliminating the imaginary portion of that sample prior to applying the term. 
In that case, the phase shift term is conjugate symmetric, and the resulting signal should be 
shifted in the time domain without incurring an imaginary portion. Results of that process 
are provided in Figure A.7, which again depicts errors relative to a signal delayed via 
interpolation. The error in the signal energy is presented in Table A.4, which shows the 
energy in the imaginary portion has been reduced to zero as expected. 
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Real part of time domain signals 



Figure A.7. Signals and error from fractional phase shifts with forced symmetry. 

Table A.4. Energy summary for fractional phase shifts with forced symmetry. 

Energy (J) 

Error 7.2528 x 1(T* 

Real(Error) 7.2528 x 10" 8 

Imag(Error) 0 

A.3 Computation Time Comparison 

If speed is a consideration, and the application can tolerate the anticipated error, delay 
by the fractional phase shift method, including a correction to eliminate the imaginary 
portion of the result, should be considered. The interpolation method requires repeated 
computation of the sinc() function, which, for a large number of samples in the signal or 
desired points of interpolation, can be computationally intensive. The phase shift method 
requires computation of the FFT and IFFT. A comparison of the times required for the 
above evaluation is provided in Table A.5. While this is just a small data set for a single 
signal, it is clear that the interpolation requires considerably more time to compute than the 
delays via phase shift in the frequency domain. Fittle more should be extracted from this 
data, given the minimal number of runs to evaluate performance. 


A.4 Summary 

When shifting an integer number of time-domain samples, a circular shift operation (e.g. 
using circshiftO in MATFAB) is a straightforward option, as this simply involves ad¬ 
justing the index numbers of the sampled data and should yield minimal, if not zero, error. 
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Table A.5. Calculation times for various shifting methods. 

Calculation time (ms) 


Circular Shift 

3.3 

Interpolation (integer) 

355.2 

Interpolation (fractional) 

358.2 

Phase Shift (integer) 

1.7 

Phase Shift (fractional) 

2.0 

Phase Shift (fractional w/ correction) 

4.4 


For fractional shifts, options include either interpolation at the delayed time or use of a 
fractional phase shift in the frequency domain. As discussed above, there is some error 
between those two methods, and while a true “control” for that error is not available, it is 
reasonable to attribute the discrepancy to the fractional phase shift method, as it involves 
a departure from the typical case while interpolation accuracy should not be affected by 
the times at which the signal is interpolated. Depending upon the application, computation 
time might be a considerable argument for use of the delay via phase shift in the frequency 
domain. 
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APPENDIX B: 

Additional Analysis Tables and Figures 


B.l OFDM Signal Recovery 


- 

.■lit 

| -0 Transmitted 

«b 

Hi u nil 

Recovered | 

O 

A h r 

O 

9 

I 

111 

0 10 

20 30 

trsec 

40 50 

S 

60 

- w©* 

o 

Its 

h i Tl 

o 

> 

fo wf 

o_? 

n 

liilllil 

^©- 1 - 


0 10 20 30 40 50 60 

Usees 


Figure B.l. OFDM test signal, time domain, parameter set 1. 
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Figure B.2. Frequency domain of OFDM test transmitted and recovered signals, parameter set 
1. 
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Figure B.3. Error of recovered OFDM signals, parameter set 1. 
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Figure B.4. OFDM test signal, time domain, parameter set 2. 
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Figure B.5. Frequency domain of OFDM test transmitted and recovered signals, parameter set 
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(a) Frequency domain (b) Time domain 

Figure B.6. Error of recovered OFDM signals, parameter set 2. 
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Figure B.7. OFDM test signal, time domain, parameter set 3. 
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Figure B.8. Frequency domain of OFDM test transmitted and recovered signals, parameter set 

3. 
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recovered OFDM signals, parameter set 3. 




Figure B.10. OFDM test signal, time domain, parameter set 4. 
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Figure B.ll. Frequency domain of OFDM test transmitted and recovered signals, parameter set 

4. 




(a) Frequency domain (b) Time domain 

Figure B.12. Error of recovered OFDM signals, parameter set 4. 


68 








































































































B.2 Geometry Analysis 


Table B.l. Recovery performance figures for layout 1. 



Signal 1 

Signal 2 

Power Angle (°) 

39.17 

48.65 

Error Figure 

0.002572— 9.5719 x 10 -29 / 

0.001851+9.2041 x 10~ 29 z 

Subspace Angle (°) 

2.907 + 3.0395 x 10 -14 z 

2.466 + 3.3244 x 10~ 13 z 

Correlation 

0.9987 - 2.6904 x 10“ 17 / 

0.9991 — 2.496 x 10 16 / 

Table B.2. Recovery performance figi 

ures for layout 2. 


Signal 1 

Signal 2 

Power Angle (°) 

39.11 

48.65 

Error Figure 

0.001286 — 3.2622 x 10~ 28 z 

0.0009253 + 1.4881 x 10~ 28 / 

Subspace Angle (°) 

2.055— 3.1089 x 10~ 13 z 

1.743 + 1.4709 x 10~ 12 z 

Correlation 

0.9994 +1.9458 x 10~ 16 / 

0.9995 - 7.809 x 10^ 16 / 

Table B.3. Recovery performance figi 

ures for layout 3. 


Signal 1 

Signal 2 

Power Angle (°) 

39.11 

48.65 

Error Figure 

0.001286 — 2.6708 x 10~ 28 z 

0.0009253 + 6.2858 x 10~ 28 z 

Subspace Angle (°) 

2.055 + 1.1063 x 10~ 13 z 

1.743 + 9.6102 x 10^ 13 / 

Correlation 

0.9994 — 6.9239 x lO -17 / 

0.9995-5.1021 x 10 - {6 i 

Table B.4. Recovery performance figi 

ures for layout 4. 


Signal 1 

Signal 2 

Power Angle (°) 

39.11 

48.65 

Error Figure 

0.001286 — 2.9969 x 10 _30 z 

0.0009253 + 2.2486 x 10 -29 z 

Subspace Angle (°) 

2.055+ 4.6904 x 10~ 14 z 

1.743 — 2.3755 x 10~ 13 z 

Correlation 

0.9994 —2.9356 x 10~ 17 z 

0.9995 + 1.2612 x 10^ 16 / 
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Table B.5. Recovery performance figures for layout 5. 


Power Angle (°) 
Error Figure 
Subspace Angle (°) 
Correlation 


Signal 1 Signal 2 


39.17 

0.002572 + 5.2147 x 10~ 23 f 
2.907 + 2.5224 x K^ 9 ^ 
0.9987 - 2.2326 x 10~ 12 f 


48.65 

0.001851 -3.0583 x 10- 21 / 
2.466-8.7371 x 10~ 10 / 
0.9991 + 6.5599 x 10~ 13 z‘ 



Frequency (kHz) 



Frequency (Mhz) 

(a) Frequency-domain error (b) Condition number 

Figure B.13. Geometry induced recovery errors, layout 1. 
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(a) Frequency-domain error 


(b) Condition number 


Figure B.14. Geometry induced recovery errors, layout 2. 
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(a) Frequency-domain error 

Figure B.15. Geometry induced 


(b) Condition number 
recovery errors, layout 3. 



(a) Frequency-domain error 

Figure B.16. Geometry induced 


(b) Condition number 
recovery errors, layout 4. 
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(a) Frequency-domain error 

Figure B.17. Geometry induced 


(b) Condition number 
recovery errors, layout 5. 
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(a) First signal (b) Second signal 

Figure B.18. Frequency domain of recovered signals, layout 1. 




(a) First signal (b) Second signal 

Figure B.19. Frequency domain of recovered signals, layout 2. 
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(a) First signal (b) Second signal 

Figure B.20. Frequency domain of recovered signals, layout 3. 




(a) First signal (b) Second signal 

Figure B.21. Frequency domain of recovered signals, layout 4. 




Frequency (kHz) 



(a) First signal (b) Second signal 

Figure B.22. Frequency domain of recovered signals, layout 5. 
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Figure B.23. Time-domain signals and error, layout 1. 
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Figure B.24. Time-domain signals and error, layout 2. 
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Figure B.25. Time-domain signals and error, layout 3. 




1 

Transmitted-© Recovered | 


- 

\,T 

J 

u 

) 

o© 

II, 

O < 

o o 

Jll 

J 

offf 
?? 1 

nl 

_ O© ffi 

iU 1 

0 

10 

2C 

30 40 

(isecs 

50 60 


1 

O 

1 

< 

< 

< 

as 

fl? 

% 

T ff ■ 
0 ° 

'll 

„ 

K 

i 

r. it 

111 Hi, 

ill, 


0 10 20 30 40 50 60 

Usees 



Time (ms) 
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Figure B.26. Time-domain signals and error, layout 4. 




(a) Time-domain signals (b) Time-domain error 

Figure B.27. Time-domain signals and error, layout 5. 
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APPENDIX C: 
Code Listings 


C.l delay s_from_lay do wn.m 


function [ tau ] = delays_from_laydown( L ) 

% DELAYS_FROM_LAYDOWN Compute the delay matrix TAU from laydown L 
% Compute signal delays from a set of transmitters to a set of 
% receivers, with geometric laydown as specified in L. 

Re = 6371e3; % m 

num.recv = length(L.r.x); 
num.xmit = length(L.t.x); 

% loop to compute all elements 
for m = 1:num.recv 

for n = 1:num.xmit 

theta = (L.r.x(m) — L.t.x(n)) /Re * 180/pi; 

% Law of cosines to get from theta to d(l,l) 
dx = sqrt ( (L. r . alt (m)+Re) "2 + Re A 2 — ... 

2*(L.r.alt(m)+Re)*Re*cosd(theta) ); 
d(m,n) = dx / cos(atan2(L.r.y(m)—L.t.y(n),dx)); 

end 

end 

% Convert from distance to time, using an approximation of the speed of 
% light 

tau = d ./ 3e8; 

end 


C.2 block_deltas_from_delays.m 
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function [ D, cond_num ] = block_deltas_from_delays( tau, omega ) 
%BLOCK_DELTAS_FROM_DELAYS Build the block diagonal matrix of delays 
% [D, cond_num] = BLOCK_DELTAS_FROM_DELAYS returns the block diagonal 

% delay matrix using the propagation delays found in TAU and the 
% frequencies from OMEGA. 

o, 

o 

% The condition number matrix allows checking for degenerate cases 

[num.recv, num.xmit] = size(tau); 

N = length(omega); 

D = exp( —lj .* reshape(tau, [num.xmit*num.recv 1]) * omega); 

D = num2cell(reshape(D, [num.recv num.xmit N]), [12]); 

cond_num = squeeze(cellfun(@cond,D)) ; 

D = blkdiag(D{:}); 

end 


C.3 calc_signals_nb.m 


function [ s, r ] = calc_signals_nb( signal_spec, t, tau ) 
%CALC_SIGNALS_NB Calculate sampled narrowband xmit and recv signals 

o, 

o 

[num.recv, num.xmit] = size(tau); 

N = length(t); 
dt = t (2 ) ; 

% Only want to preserve the fundamental freqs and magnitudes 
s = struct (' mag : , signal_spec.mag, 'f', signal_spec.f); 

for n = 1:num.xmit 

s.func{n] = @(t) (s.mag(n)*cos(2*pi*s.f(n)*t)) ; 
s.time(n, :) = s.func{n}(t); 
s.freq(n,:) = fft(s.func{n}(t))*dt; 

end 
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% Received signals — "r" 
r.time = zeros(num.recv,N); 
for m = 1inum.recv 

for n = l:num.xmit 

r.time(m,:) = r.time(m,:) + s.func{n}(t — tau(m,n)); 

end 

end 

r.sig_power = sum(r.time. * 2 , 2)/N; 
end 


C.4 recover_from_delays.m 


function [ s_recovered ] = recover_from_delays( num, D, sig_freq, dt, 
error_check) 

%RECOVER_FROM_DELAYS Use signal delay to recover original signal 

o, 

o 

% Recover the original signals from a set of co—channel signals 
% recevied at a separated set of receivers, given the known delays 
% between the transmitters and receivers. 

o, 

o 

% s_recovered = recover_from_delays( num, D, sig_freq, dt, 

% error_check ) returns the original signals in s_recovered. D is a 

% mN x nN matrix of the signal delays where m is the number of 
% receivers, n is the number of transmitters, and N is the number of 
% discrete frequencies. sig_freq is the matrix of frequency 
% samples, generally the result of fft(time_domain_samples). dt is 
% the sampling interval in the time domain and error_check is a 

% boolean that controls an extra level of error checking in the 

% recovery process, useful for understanding errors at discrete 
% frequencies due to poorly conditioned blocks in the D matrix. 

N = length(sig_freq); 

R = reshape(sig_freq, [num.recv*N,1]); 

% remove the rows and columns that correspond to omega(1) = 0 
D = D(num.recv+1:end,num.xmit + 1:end) ; 
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R = R(num.recv+1:end); 


% preallocate S for speed and proper dimension, 

% necessary to account for the missing elements per above 
S = zeros(num.xmit*(N-l), 1) ; 

if error_check 

% cycle through the blocks to ensure independent calculation 
% there are N—1 blocks, after the removal completed above 
for blk_num = 1:N—1 

index = 2*(blk_num—1) + (1:num.xmit); 
if cond(D(index,index)) > 300 

fprintf('Bad condition number at indices: [%d %d]\n', ... 

index+num.xmit) %adjust for missing block 
disp(D(index,index)) 

end 

S(index) = D(index,index) \ R(index); 

end 

if (-all(D * S — R < eps) | -all(D \ R — S < eps)) 

warning( 1 Recovery had errors greater than epsilon.') 

end 

else 

S = D \ R; 

end 

S = reshape(S, num.xmit, []); 

% Insert zeros for the known singularity at f=0 and return to time domair 
s_recovered.freq = [zeros(num.xmit,1) S] ; 
s_recovered.time = ifft(s_recovered.freq, [], 2)/dt; 

s_recovered.sig_power = diag(s_recovered.time * s_recovered.time') / N; 

s_recovered.N = length(s_recovered.time); 

end 


C.5 calc_signals_wb.m 


function [ s, r ] = calc_signals_wb( signal_spec, tau ) 
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CALC_SIGNALS_WB Calculate sampled wideband transmit and receive 
signals 

[ s, r ] = CALC_SIGNALS_WB(SIGNAL_SPEC, TAU) returns the 
transmitted (S) and received (R) signals, with delays specified by 
TAU incorporated in the received signals. 

SIGNAL_SPEC is a structure that includes SIGNAL_SPEC.TYPE, which 
may be either 'noise' or 'impulse' for a AWGN signal, or a time 
domain impulse, which has a constant response in the frequency 
domain. The impulse is useful when testing. SIGNAL_SPEC.N is the 
number of points, and SIGNAL_SPEC.dt is the inter—sample interval, 
or 1/Fs. 

CALC_SIGNALS_WB() with SIGNAL_SPEC.TYPE = 'noise' returns a signal 
with power spectral density equal to 1. 

If defined, SIGNAL_SPEC.RNG_STATE is used to set the state of the 
RNG prior to signal generation. This ensures consistent signal 
contents across multiple calls to CALC_SIGNALS_WB(). 


s = signal_spec; 
dt = signal_spec.dt ; 

% Use a seed for the RNG, if requested 
% constant 

if isfield(signal_spec, 'rng_state') 
rng(signal_spec.rng_state) 

end 


keeps signal contents 


[num.recv, num.xmit] = size (tau); 

buffer_len = ceil(max(reshape(tau, [num.recv*num.xmit 1])) / dt); 


for n = 1:num.xmit 

switch (lower(signal_spec.type)) 
case 'noise' 

N = signal_spec.N; 

N_hat = 2. A nextpow2(N+buffer_len); 
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s.time(n,:) = randn(l,N); 

s.freq(n,:) = fft(s.time (n, : ) , N_hat, 2)*dt; 

case 'impulse' 

N = signal_spec.N; 

N_hat = 2. A nextpow2(N+buffer_len); 
s.freq(n,:) = ones(1,N_hat); 

s.time(n,:) = ifft(s.freq(n ,:), N_hat, 2)/dt; 
case 'ofdm' 

signal = ofdm_signal(s.ofdmParams{:}); 

% The OFDM signal needs to be interpolated dt time domain 
N = ceil(signal.duration / dt); % new number of samples 

N_hat = 2. A nextpow2(N+buffer_len); 

t = (0:N_hat—1)*dt; % new time domain, based on dt 

s.time(n,:) = interp_time(t, signal.sig, signal.Tsamp); 
s.freq(n,:) = fft(s.time (n, : ) , N_hat, 2)*dt; 

end 

end 

df_hat = 1 / (N_hat*dt); 

fprintf( 's.time energy: %f\t%f\n', diag(s.time * s.time')*dt) 
fprintf( 's.freq energy: %f\t%f\n' / diag(s.freq * s.freq')*df_hat) 

o, o. 
o o 

% create a new omega vector for the longer time vector 

k = (0:N_hat—1) — N_hat/2; 

omega = fftshift(2*pi*k / (N_hat*dt)); 

o, o. 
o o 

% Delayed signals (as though received independently at the receivers) 
for m = linum.recv 

for n = l:num.xmit 

delayed.freq(m,n,:) = s.freq(n,:) .* exp(—1j*omega*tau(m,n)); 

end 

end 

delayed.time = ifft(delayed.freq, N_hat, 3)/dt; 
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for 


m 


1:num.recv 


for n = linum.xmit 

del.t = squeeze(delayed.time(m,n,:)); 
del.f = squeeze(delayed.freq(m,n,:)); 
fprintf (' delayed.time(%d,%d) energy: %g\n', ... 

m,n,diag(del.t' * del.t)*dt) 
fprintf (' delayed.freq(%d,%d) energy: %g\n*, ... 

m,n,diag(del.f' * del.f)*df_hat) 

end 

end 

o, o. 
o o 

for m = 1:num.recv 

for n = linum.xmit 

corr_plot(s.time(n,:), ... 

permute(squeeze(delayed.time(m,n, :)) , [2 1])); 

end 

end 

o, o_ 
o o 

% Received signals — "r" 
r.time = zeros(num.recv,N_hat); 
for m = 1:num.recv 

for n = linum.xmit 

r.time(m,:) = r.time(m,i) + ... 

permute(squeeze(delayed.time(m,n,:)), [2 1]); 

end 

end 

%r.time = r.time(:,1iN); 

r.freq = fft(r.time, N_hat, 2)*dt; 

fprintf(' r.time energy: %f\t%f\n', diag(r.time * r.time')*dt) 
fprintf(' r.freq energy: %f\t%f\n', diag(r.freq * r.freq')*df_hat) 


o, o_ 
o o 

% Signal powers 

s.sig_power = diag(s.time * s.time') /N_hat; 
r.sig_power = diag(r.time * r.time') /N; 
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end 


C.6 interp_time.m 


function output = interp_time( t, x, Tsamp ) 

% INTERP_TIME Interpolate using the ideal low pass filter. 

% INTERP_TIME( tO, x, Tsamp ) returns an interpolated estimate of 
% x(tO) where tO is the time of the desired point for interpolation, 
% x contains samples of the function x() such that x[n] = x(n*Tsamp) 
% and Tsamp is the sampling interval. 

o, 

o 

% INTERP_TIME( t, x, Tsamp ) returns a column vector of interpolated 
% values, where interpolation is performed at times specified in t. 

o, 

o 

% Authored by Ned Patterson 

n = 0 : length(x)—1; 
n = ones(length(t),1) * n; 

t = reshape(t, [], 1); % Ensure that t is a column vector 

t = t * ones (1, length(x)); 

% Assume we have sampled at the nyquist rate 
B = 1/(2*Tsamp); 

h = sine( t/Tsamp — n ); 

x = reshape(x, [], 1); % Ensure that x is a column 

% Compute the sum of the product using the inner product 
output = h * x; 
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